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ABSTRACT 

> 

\^ _ Aims. We devise a data reduction and calibration system for producing highly-accurate 21 -cm 

; HI spectra from the Green Bank Telescope (GBT) of the NRAO. 

Methods. A theoretical analysis of the all-sky response of the GBT at 21 cm is made, augmented 
by extensive maps of the far sidelobes. Observations of radio sources and the Moon are made to 
check the resulting aperture and main beam efficiencies. 

Results. The all-sky model made for the response of the GBT at 21 cm is used to correct for 
"stray" 21-cm radiation reaching the receiver through the sidelobes rather than the main beam. 
This reduces systematic errors in 21-cm measurements by about an order of magnitude, allowing 
' accurate 21-cm HI spectra to be made at about 9' angular resolution with the GBT. At this 

resolution the procedures discussed here allow for measurement of total integrated Galactic H I 
line emission, W, with errors of 3 K kms"', equivalent to errors in optically thin A'hi of 5 x 
10'^ cm-\ 

Key words. Radio lines: ISM - Methods: observational - Methods: data analysis - 
Instrumentation: detectors 

1. Introduction 

Accurate spectra of Galac tic HI are needed for many areas of research , from studies of inter- 



stellar gas and dust (e.g., iHeiles et alJ Il981t iBoulanger & Peraulti Il988c iBoulanger et al. 



1996 
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Lockman & Condo 
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Kalberla & Kert]l2009t 



Planck Collaboration et al 



nation of the distribution and abundances in interstellar gas (e.g. 



Hobbs et al 



201 lah . to determi- 



1982; 



Albert et al 



1993 



ShuUetal 



20091) . to corr ection of exlragalactic radiation for absorption by the ISM 



and removal of foregrounds (e.g. , 



Hauser et al 



1998 



Puccetti et al 



Jahoda et al 



2011 



1985 



Hasinger et al 



Planck Collaboration et al 



1993 



Snowden et al 



1994 



201 Ibh . Although the highest 



resolution H I observations are obtained with aperture synthesis techniques, Galactic H I emission 
is smoothly distributed across the sky, i.e., most of the power is in the lowest spatial frequen- 
cies, and so filled aperture (i.e., single dish) observations are essential for determining accura te H I 
spectra (e.g.. 



Heiles&Wrixonll976: 



nll976[lGreeJl993l 



Dickev et alj|2001 



Kalberla & Kerdl2009 ). The 



i 



most accurate single-dish Galactic H I measurements are now capable of determining W, the inte- 
grated emission over the line profile, with an error of just a few percent, and hence yield A^hi, the 
total column density of Galactic H I to the same precision provided that opacity effects are small 
Jwakkeret allboilL 

The Robert C. Byrd Green Bank Telescope (GBT) is a 100-m diameter fiUed-aperture radio 
telescope that has been used for many studies of Galactic H I and the extended H I emission 
around nearby galaxies (e.g.. iHunter et al.ll201 ih . In this paper we describe techniques whereby 



we are able to produce high-quality 21 -cm spectra with this instrument with overall calibration 
errors of only a few percent. For illustration, we consider data from rece nt HI surveys under- 



taken to study the gas-dust correlation and dynamics of high latitude cirrus ( IBlagrave et al 



2O10; 



Planck Collaboration et al.ll201 la : 



Martin et al.ll201 11) . although the methods described here can be 



applied to any 2 1 -cm H I data taken with the GBT. It is important to note that measurement of accu- 
rate extragalactic H I is much more straightforward than for Galactic H I as there is little confusing 
emission entering through sidelobes at |v| > 200 km s ' ; errors in extr agalactic H I profi les can be 
as small as 3%, arising largely from instrumental baseline effects (e.g., Hogg et aklboO? ). 

In Sect. |2]we describe the telescope and its instrumentation. Section[3]considers the observing 
techniques used for our Galactic H I measurements. In Sect. |4] we describe the steps taken to reduce 
and calibrate the data. Section |5] describes a theoretical calculation of the main beam properties, 
the caUbration of the antenna temperature scale, conversion to brightness temperature, and several 
independent cross-checks of the accuracy of these steps. In Sect. |6] measurements of the antenna 
response up to 60° from the main beam are described, leading to an all-sky antenna response 
suitable for correcting for "stray radiation" detected through the sidelobes rather than the main 
beam. Section l63] illustrates the effects of the correction on several 21 cm spectra. In Sect.|2]several 
tests of the data reduction method and our estimates of the errors are documented. In Sect. [8] we 
discuss the absolute calibration. Section|9]presents a summary and conclusions. 



2. The Green Bank Telescope and its instrumentation 



The GBT is a 100-m diameter dual offset Gregoria n reflector with a larg e, unblocked aperture on 



an azimuth-elevation mount (Fig. [T] Fig. |2l Fig. [3j IPrestage et al 



20091) . The surface consists of 



2004 panels mounted on motor-driven actuators capable of real-time adjustment to compensate for 
gravitational astigmatism and other surface distortions allowing it to achieve a surface rms < 250//. 
We used the telescope with the surface in passive mode where it has a typical rms error =s 900/v, and 
thus at 21 cm (1420 MHz) a typical surface efficiency of 0.997. Because of gravitational distortions 
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Fig. 1. The Green Bank Telescope in a view that shows its unblocked 100-meter diameter aperture. 



and thermal effects, the passive surface rms accuracy of the GBT is expected to vary between SOOyu 
and 1200/1 at the extremes. This would cause a slight change in the main beam shape and a point- 
source gain change with a range of up to 0.4%. The effect on HI observations is more complex 
and likely less important as it depends on the convolution of the detailed beam shape with the H I 
sky brightness. The main beam has a FWHM that is 9.1' x 9.0' in the cross-elevation and elevation 
directions, respectively. At 21 cm, the aperture efficiency rj^ = 0.65 and the main beam efficiency 
TJmb - 0.88. The derivation of these quantities will be discussed in Sect. 15. H and AppendixlAl 

The 21 -cm "L-band" receiver used for these measurements is located at the secondary 
Gregorian focus and illu minates an 8-m diameter subreflector with an average edge taper of 
-14.7 dB jSrikanthlll993l) . The receiver is cryogenically cooled, accepts dual linear polarization, 
and has a total system temperature, Tsys, at the zenith of 18 K. A diode injects a fixed amount of 
noise into the waveguide just after the feed horn, before the polarizer and amplifier, and can be 
modulated rapidly for calibration. This was used to establish a preliminary scale for the antenna 
temperature, T^. These noise sources are typically very stable; we see no evidence for variation of 
the L-band receiver calibration diode over several years (Sect.|7]i. 



3. H I 21 -cm line data acquisition 

3.1. Mapping 

The H 1 2 1 -cm observations discussed here were made with On-the-Fly (OTP) mapping: "scanning" 
or moving the telescope in one direction, typically Galactic longitude or Right Ascension, while 
taking data continuously. The integration time and telescope scan rate must be chosen so that 
samples are taken no more coarsely than at the Nyquist interval, ^ FWHM/2.4 = 3.8' for the 
GBT at 21 cm, and ideally at half that interval to avoid beam broadening in the scanning direction 
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Fig. 2. The focal area of the GBT as seen from the surface of the main reflector, showing outUned 
in red, bottom to top, (i) the receiver room, with (ii) the L-band feed horn pointing upward at (iii) 
the 8-meter subreflector. Just below the lower edge of the subreflector the lower portion of a screen 
is visible. This functions to direct feed spillover radiation that would otherwise strike the telescope 
arm back into the main reflector and onto the sky. 



dMangum et alj2007h . Areas were mapped by stepping the scans in the fixed coordina te (the cross 



scan direction) and reversing the scan direction. In the HI surveys discussed here (iMartin et al 



201 ll) scans up to 5° long were made, but for practical purposes larger regions were broken up into 
smaller areas with dimensions between 2° x 2° and 4° x 4° that were mapped independently. 



3.2. Spectrometer setup 

The GBT autocorrelation spectrometer as configured for these measurements has 16k channels over 
a 12.5 MHz band with nine-level sampling in each linear polarization, for a velocity resolution of 
0.16 kms ' in the 21 -cm line. Early in the data reduction procedure the resolution was reduced 
to 0.80 kms ' by filtering the spectra with an eleven-channel Hanning smoothing function and 
resampling every fifth channel. This maintains the independence of each channel, while minimizing 
aliasing effects. 

The local oscillator was modulated to move the center of the spectrometer band between a 
"signal" and "reference" frequency, separated by 2.5 MHz (528 kms '). The noise source was 
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Fig. 3. The GBT subreflector and feed arm seen from the secondary focal point atop the receiver 
room. The screen redirects feed spillover away from the arm and down into the main reflector. 

modulated synchronously with the frequency switching to calibrate the receiver gain at both signal 
and reference frequencies every second. Because the modulated separation is much smaller than the 
12.5 MHz covered instantaneously by the spectrometer, the emission in the 21-cm line was always 
being observed. This "in-band" frequency switching gives a factor of two increase in observing 
speed over "out-of-band" frequency switching, and an rms noise in antenna temperature, o-j.^, of 
approximately 0.25 K in a 1 kms ' channel in 1 second for the average of the two polarizations. 
For the basic integration time for the survey data, 4 seconds, and velocity resolution 0.80 kms \ 
the resulting rms noise is 0.16 K for emission-free channels (see Sect j7.1l ). Some regions were 
measured several times to check the reproducibility and to improve the sensitivity. 

3.3. Spectral baselines 

The frequency-switched spectra produced by the GBT spectrometer (e.g.. Fig. lU have remarkably 
flat instrumental baselines over the central 4 MHz, even more so for the XX compared to the YY 
polarization. Nevertheless, an important step in our GBT data reduction procedure (Sect. |4]i is 
removal of any residual instrumental baseline. As is usual, the baseline is approximated by a low- 
order polynomial whose coeflicients are fit by linear least-squares to the "emission-free" channels 
of each spectrum. For the GBT spectra we found that over the central ~ 4 MHz (-450 < v < H-400 
km s"') a third-order polynomial was adequate, the next order term being not statistically justified. 
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Fig. 4. Brightness temperature spectrum from the North EcUptic Pole (NEP) data cube at (/, b) - 
(90^3 18, 34'?432), illustrating aspects of the baseline removal. The negative feature at 370 kms"' 
results from the 2.5 MHz frequency-switching: it is an offset half-amplitude inversion of the 
- 170 km s~' feature. The upper panel shows a typical third-order polynomial baseline fit to the cen- 
tral 4 MHz of the frequency-switched spectrum; the horizontal Une segments mark the "emission- 
free" channels determined iteratively to be part of the baseline (see text). The lower panel displays 
the spectrum following baseline removal. 



Except for diagnostic purposes, we fit baselines to spectra after correction for stray radiation and 
assembly into a data cube (Sect. I4.2l i. In a very few spectra, out-of-band interference results in 
unsalvageable spectra with high-order polynomial baselines. These are easily identified by eye, and 
in most cases we were able to reobserve these positions. Alternatively, these spectra are flagged for 
exclusion from the map-making process. 

For our mapping observations with the GBT, the instrumental baseline was found to change 
slowly with time, so that the coefficients of the polynomial fit are highly correlated between spectra 
along a scan. For noi sy spectra there are p otential advantages to fitting a baseline to the average 
of sequential spectra dLockman et al.lll986h . but we did not implement this as fitting of individual 
polynomials gave adequate results. It is possible to take advantage of an iterative baseline fitting 
technique. This brings in more channels than one normally obtains from a conservative estimate 
of the location of the emission-free end channels for a particular map, even locating emission-free 
parts of the spectrum between emission features. 

We fit baselines following the iterative technique used by 



Hartmann 



(119941) for the 



Leiden/Dwingeloo Survey. Prior to any fitting, each spectrum is smoothed by a 20-channel (roughly 
16 km s ') boxcar to accentuate real velocity features. We used Hartmann's definition of a "veloc- 
ity feature" present in a residual spectrum from which an estimated baseline has been removed. 
These are found by first identifying a significant (4. Oct) positive peak and then adding neighbouring 
channels until the residual goes negative; ten additional channels (about 8 kms ') at both ends of 
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each velocity feature are also flagged for omission. Following the identification of all such features, 
there remain the "emission-free" baseline channels to be used in the next fit. The smoothed baseline 
spectrum is fit successively with a series of polynomials monotonicaUy increasing from linear to 
third-order, after every iteration subtracting this updated baseUne and identifying and flagging new 
velocity features with significant residual peaks. Because the spectra are frequency-switched, any 
velocity emission feature will show up as a half-amplitude inverted feature ofi'set by the switching 
frequency, 2.5 MHz. These channels are also flagged for exclusion. Finally, a third-order polyno- 
mial is fit to the remaining list of emission-free channels in the original spectrum and subtracted. 
Usually ~ 600 emission-free channels are used for the fit. Figure|4]illustrates a typical result of the 
iterative baseline fitting process, for a typical instrumental baseline. 

Occasionally we found H I emission from background galaxies in our spectra, usually in the end 
channels beyond the Galactic emission, but sometimes even overlapping it. This makes determining 
the instrumental baseline challenging. While this can be treated on a case-by-case basis, these pixels 
are simply masked as unsuitable for analysis of Galactic H I. 

3.4. Treatment of radio frequency interference (RFI) 

The main source of RFI in these data is an oscillator in the GBT receiver room which can produce 
narrow-band spurious signals in the data. These are extremely stable and << 1 kHz in width. They 
can be identified easily and removed. A series of eight frequency ranges in which RFI had been seen 
were systematically inspected by averaging many spectra together in the topocentric velocity frame 
of reference. Any 3.5cr positive deviations from the median over the suspect frequency range were 
flagged. Finally data in a group of five channels around the flagged frequency were replaced with 
values from a linear interpolation of the surrounding channels. This channel replacement is done at 
the highest velocity resolution (0.16 kms"') on the four observed spectrometer phases (sig/calon, 
sig/calofi", ref/calon, ref/caloff) prior to producing the final calibrated frequency-switched spectra 
and prior to any subsequent spectral smoothing. Efi'orts are underway to replace the interfering 
device. 

4. Data reduction 

The measured quantity, T^, is related to the 21 -cm H I sky brightness by 



where e is the direction-dependent atmospheric extinction and P is the antenna power pattern; 
the integration is over the entire 4jt sr and this "all-sky" integral of P is unity (later we also refer 
to the antenna gain G = 4nP relative to isotropic, where by definition Gisotropic(0, 0) = !)■ The 
quantity rji- accounts for resistive losses in the system which are less than 1%. Our procedure in 
Sect. l5.2l for calibration of Ta produces an antenna temperature scale that is independent of t/i and 
we will thus drop this factor from subsequent equations. Equation [T] can be separated into terms 
that come from the main beam and from elsewhere on the sky: 




(1) 




(2) 
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where e"'""* is the atmospheric extinction in the direction of the main beam, T/mb, the main beam 
efficiency, is the fraction of the total power accounted for by the main beam over Qmb, f^si is the 
area of the sky outside the main beam (Osi = 4;r - Q.mb)- The desired quantity to be measured, 
(T^mb), is the H I brightness temperature averaged over the main beam: 



<r„,b) 



r rb(0,<^)dQ„b. 



(3) 



While Tb and are functions of frequency because of Doppler shift, it is assumed that the 
other quantities are constant over the relatively narrow frequency range of Galactic HI emission. 
Equation |2] can thus be written 



<r„b> 



'7mb 



/ 



P(e, 0)rb(6', 0)e"^"dQ 



(4) 



The aperture efficiency 77a and beam efficiency zy^b are not necessarily independent (e.g.. 



Goldsmith 



20021). For a dish of diameter, D, at wavelength. A, with a Gaussian main beam. 



%WHM, and with typical values for the main reflector edge taper, we have T]„h = Omb/f^a, 
□a = 4A^/(i],nD^), ilmb ^ 1-130^HM' and 0fwhm * 1.25^/D. Thus, 77„,b * 1.477a for the GBT at 
21 cm. However, this relationship is not necessarily precise for real systems, and so we treat the 
two terms as independent quantities that have to be derived and checked separately. 

For simplicity, the observable (Tmb) will be referred to as T^b or simply T from henceforth. 



4.1. Stray radiation 

The second term in Eq. |4] accounts for H I emission that enters the receiver through the sidelobes 
of the telescope rather than through the main beam. This is called "stray" radiation. Just as the 
main beam efficiency, T/^b, reflects the fraction of the power pattern in the main beam, we define a 
sidelobe efficiency, rj^i, 



77si= r m0)dO 



1 



(5) 



The unblocked aperture of the GBT eliminates the scattering sidelobes that plague most other 
radio telescopes, but at 21 cm the relati vely low edge taper of the feed (originating from mechanical 
limitations on its size and weight, see 



Norrod & Srikanth 



19961) results in a "spillover lobe" past 



the secondary reflect or. Because there is 21-cm H I emission fr om all directions of the sky at a level 
A^Hi ^ 4x 10'^ cm"'^ ( Lockman et alj 1986[|jahodaetalJl99ol) . there will be a contribution to every 



spectrum from H I emission entering through the sidelobes, and so the "stray" radiation spectrum 
must be calculated and removed as a part of data reduction and calibration. Prior to construction 
of the GBT the only unblocked radio telescope useful for 21-cm measurements was the Bell Labs 



horn-reflector which had a m ain beam size ^ 3° x2° and a spillover lobe as well jWrixon & Heiles 
Kuntz&Danlvlll992i 



1972t 



Stray radiation originating from large angles is Doppler-shifted with respect to the direction of 
the main beam, so that the stray radiation spectrum correctly shifted to the local standard of rest 
(LSR) depends not only on the locati on of the sidelobe on the sky but also on the date and time of 



the observation jKalberla et al 



20051) . For an Alt-Az telescope like the GBT, the sidelobe pattern. 



which is fixed in telescope coordinates, also rotates on the sky about the main beam as the Local 
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Sidereal Time (LST) changes. A direction observed at varying LST will have a varying component 
of stray radiation. 



St ray radiation has been a n issue for Galactic 21 -cm science for 50 years dvan Woerden et al 



1962; 



(e.g.. 



Heiles & Wrixon 



Wrixon & Heiles 



Higgs et al. 



2005 



1976 ) and many groups have devised methods to 



1972; 



Kalberla et al 



Kalberla et al 



198C; 



Lockman et al. 



1986 ; 



suppress or rernove it 



Hartmann et al 



1996 



2010t) . The most accurate methods solve Eq. |4]for (r,j,b) using 



models of the antenna power pattern, P, and all-sky T\,(6, (p) maps. 

Because the most important sidelobes of the GBT cover large areas of the sky, relatively low an- 
gular resolution 2 1 -cm surveys can be used for t he Ti,(6, (p) terrn in es timating the stray component, 
e.g., the Leiden/ Argentine/Bonn (LAB) survey jKalberla et alJl2005L However, this is not true for 
the first few diffraction sidelobes that lie close to the main beam. These have angular structure on 
scales like that of the main beam and would require knowledge of the H I sky at that level of detail 
for their removal. For the GBT, however, the near sidelobes are quite small and contain < 1 % of 
the telescope response (Sect. lSTI ). We do not correct for these, in effect assuming that their small 
component of the telescope's response samples nearly the same emission as the main beam. 



4.2. Procedure 

The data reduction procedure thus involves calibration of the intensity scale to antenna tempera- 
ture, calculation and subtraction of the stray radiation spectrum, and correction for the main beam 
efficiency and the atmosphere. Data reduction additionally requires removal of an instrumental 
baseline, and for maps, interpolation of the sampled spectra into a data cube. The latter two steps 
are sometimes interchanged; the interpolated spectra, especially if combining several observations 
of a region, are less noisy so that baseline removal at that final stage might be preferable. 

The initial calibration of the data to an approximate antenna temperature scale used values 
for the receiver calibration noise source that were determined by measurements in the laboratory 
for each receiver polarization channel. This part of the data reduction was performed using the 
NRAO program GBTIDL. A constant calibration temperature was assumed over the 12.5-MHz 
band. We checked this part of the calibration, and we correct it by a small amount using our own 
measurements and calculations as described below in Sect. 15.21 

Our implementation of the stray radiation correction is described in Sect. 16.41 and AppendixICl 
making use of the sidelobe pattern established in Sect. 16.11 

Calculation of the main beam efficiency is described in Sect. 15.11 The NRAO data reduction 
program GBTIDL will perform an approximate correction for atmospheric extinctioij^but we chose 
to make the correcti on independently , usi ng Tzenith = 0.01036 + 0.00059, the weighted mean of the 
measured values of IWilliamsl (Il973h and Ivan Zee et al.l (Il997h . with a model for the atmospheric 



air mass (AppendixO. 

After the spectra were calibrated and corrected for stray radiation, an instrumental baseline 
was removed from each spectrum by fitting a third-order polynomial to emission-free velocities. 
It is important that the stray radiation be removed before baseline fitting, lest weak stray wings 
be mistaken for instrumental baseline. For mapped regions, data cubes were constructed in classic 



' Jim Braatz, October 30, 2009: Calibration of GBT Spectral Line Data in GBTIDL v2.1, from 
|http://www.gb.nrao.edu/GBT/DA/gbtidl/gbtidl-calibration.pdf| 
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Table 1. Calculated forward gain and aperture efficiency for the GBT L-band receiver 



Frequency 


Forward Gain 


'7a 


(GHz) 


(dBi) 




1.15 


60.1 


0.708 


1.40 


61.6 


0.654 


1.73 


63.7 


0.712 



AIPS, averaging togetlier spectra from the two polarizations, using the optimal tapered Bessel 
function for interpolation (IMan^um et alJl2007h . In our data reduction procedure, baselines were 
removed after creating the cubes. For details see Sect. 13.31 



Note that nowhere du ring this procedure are "standard" H I regions, like S6 (IWilliamsll973h and 
S8 ( Kalberla et al.l[l982 ). used to calibrate the GBT intensity scale. Because of the clean optics of 
the GBT, we preferred to determine its characteristics from basic calculations and radio continuum 
flux density calibration sources. We did, however, observe the above two standard H I directions 
and we discuss these measurements in Sects. 1741 and l8. II We also compared our spectra to those 
of the LAB survey (Sect. [ 



5. Aperture and main beam efficiencies 

5.1. Calculation of the efficiencies 

A theoretical estimate of the all-sky response of the GBT was calculated using a reflector antenna 
code developed at the Ohio State University as described in Appendix lAl This incorporates infor- 
mation on the detailed illuminatio n pattern of the 21 -cm feed on the GBT subreflector as measured 
after construction of the receiver JsrikanflJll993L Table[T]shows the calculated on-axis or forward 
gain and aperture efficiency as a function of frequency for the L-band receiver. The variation in 
results from measured changes in the L-band receiver illumination pattern with frequency. The 
calculated gain G{0, (p) of the GBT at 1 .4 GHz within 1°2 of the main beam is shown in Fig.|5]for ra- 
dial cuts along polar angle in planes at several angles (p (see the angle definitions in AppendixlAb. 
The units are dBi (logarithmic units relative to an isotropic beam) and the value at = 0° is that 
in Table [1] The first sidelobe is calculated to be 29 dB below the forward peak of the main beam. 
Observations also indicate that the GBT's main beam is exceptionally clean with near-in sidelobes 
all about 30 dB below the main beam gain jRobishaw & Heilesll2009l) . 

Table |2] gives the fraction of the total power pattern lying within a given radius around the 
main beam, at radii corresponding to the minima in Fig.|5] This shows that 87.7% of the antenna 
response lies within Q°.4 of the main beam, increasing by only a small amount to 88.1% at r.2 
radius, consistent with the observations that the near sidelobes are at very low levels. We adopt a 
value 77mb = 0.88. This valu e was also derived independently during the initial calibration of the 
L-band receiver on the GBT JHeiles et alJl2"()03i . 

We define the "far" sidelobes as those arising at angles > 1° from the main beam. Compared 
to 77si 0.1, the potential for stray radiation arising within (f.4 - 1° of the main beam is negligible 
(At; » 0.004); furthermore, the brightness of this radiation T\, will not be grossly different from that 
seen on axis, and it will not be Doppler shifted as in the far sidelobes. 
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Table 2. Calculated fractional power pattern at 1 .4 GHz as a function of angle from the beam center 



e 


P(0) 


ffi26 


0.871 


a4i 


0.877 


K55 


0.880 




0.881 


0:79 


0.881 


a 89 


0.881 


1°20 


0.881 



CD 



70 



60 



50 



40 



30 



m 20 
10 





-10 



-20. 



-1.5 



-1.0 



-0.5 



0.0 

0° 




0.5 



1.0 



1.5 



Fig. 5. GBT antenna gain above isotropic for the main beam and near sidelobes at 1 .4 GHz, as a 
function of polar angle 6 from the main beam in several planes defined by cp (see Appendix lAl for 
the antenna code used for this calculation). The forward gain is 61.6 dBi (Table[T)- 



5.2. Establishing fhe antenna temperature scale 

The preliminary antenna temperature calibration assumes that the noise source has no frequency 
dependence over the range of our observations. To check this we observed the standard radio con- 
tinuum source 3C286 in spectral line mode and measured the precise value of the receiver noise 
source averaged in 50 kms ' intervals around the Galactic 21 -cm line. When the two linearly po- 
larized channels are averaged the noise source varies little with frequency, with a peak-to-peak 
fluctuation about the mean of only 1 .2% over 600 km s"' around the 21-cm line. Thus the assump- 
tion of a constant caUbration noise across our band will not contribute a significant uncertainty to 
the HI measurements when both polarizations are combined. 

The absolute value of the noise source was checked using the flux density standard 



3C286 (lOtt et al.l 119941) and the aperture efliciency derived from the electromagnetic calcula- 
tions (Sect. l5Tt . From 13 measurements at two epochs we derive a calibration temperature of 
1.495 + 0.013 K (I0-), a value that is in the ratio 1.024 + 0.009 to that measured in the laboratory. 
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As 3C286 has significant Unear polarized emission this is for the average of the two polarizations. 
A smaller number of measurements on 3C295 gives a result that is identical to within the uncer- 
tainties. We thus increased the preliminary antenna temperatures by this ratio to place spectra on 
an accurate T-^ scale. Note that this calibration, derived from an external radio source, subsumes 
within it any correction necessary for ohmic losses in the entire system. 



5.3. A check using the Moon 
The Ta scale and value of 77n,b were further checked by continuum observa tions of the Moon, 



assumed to have a constant brightness temperature of 225 ± 5 K at 1 .4 GHz (IKeihm & Langseth 



1975b and taking into account that only a part of the response pattern shown in Fig.|5]lies within the 
disk of the Moon. Measurements at several epochs give a measured to expected ratio 0.97 + 0.025 
(Icr), where the uncertainty includes the scatter in the measurements and in the assumed T\, of 
the Moon. Because we measure the Moon relative to the nearby sky, these observations are not 
susceptible to radiation in the far sidelobes and so are an uncomplicated test of both the scale 
and ?7mb determination. We conclude that our calibration does not have systematic errors that exceed 
a few percent. 



6. GBT all-sky response 

Calculating the stray 2 1 -cm component requires knowledge of the antenna response in all directions 
on the sky. Important sidelobes can occur at a level 50 dB below that of the main beam, and thus 
be quite difficult to measure. Nonetheless, they can contain several percent of the telescope's total 
response because they cover a large area on the sky . In an analysis of the 21 -cm beam pattern of 



the EfFelsberg 100-m antenna. iKalberla et al.l d 19801) found 70% of the response within 15' of the 



main beam, another 12% in the range 15' < 6 < 4°, and the remaining 18% at 4° < < 180°. 
An unblocked antenna like the GBT has much lower sidelobe levels, but at 21 cm there is still an 
important contribution to stray radiation due to spillover past the subreflector. Note that the sidelobe 
pattern for the GBT is not symmetrical about the main beam, but for comparison more than 87% of 
the response is within 15' of the main beam (Table|2]l, with another 1.8% in the range 15' < 6* < 4°. 

To determine the all-sky response P of the GBT we first used the reflector antenna code 
(Sect. 15. Il l to find the calculated response. The results were used to set the aperture efficiency, 
main beam efficiency, and the power pattern within 6 < 1° of the main beam. 

These calculations also give a general picture of the location and amplitude of the far sidelobes. 
Key features of the sidelobes can also be understood using simple near-field diffraction theory as 
discussed in Appendix |B] Because the antenna code uses only an approximation to the complex 
structure of the GBT, a more accurate determination was accomplished by measuring the sidelobes 
directly over much of the region within 60° of the main beam, exploiting the Sun as a strong radio 
continuum source. The Sun can be treated as a point source relative to the angular size of most of 
the structure in the sidelobes. 
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Fig. 6. Mapping the GBT sidelobe pattern using the Sun. Heavy hnes indicate the position of the 
Sun relative to the main beam during the scan. Thin magenta hnes indicate scans reflected in the 
beam's hne of symmetry (i.e., with the H coordinate replaced by -H). Dotted lines at high H show 
measurements that were not used because they showed no visible sidelobes above the noise. The 
center of the subreflector is 12.3° above the main beam. 



6.1. Sidelobe measurements using ttie Sun 
6.1.1. Observations 

The observations were made on three diff'erent occasions, during times of the year when confusion 
with radio emission from the Galactic plane would be minimized. Data were taken over a 20 MHz 
band centered on 1420 MHz. The mapping procedure consisted of raster scans, moving the tele- 
scope either in elevation ("vertical" scans) or azimuth ("horizontal" scans). Figure |6] shows the 
position of the Sun relative to the main beam during the scans used to determine the sidelobes, in 
terms of the vertical separation in elevation V from the beam center, and the horizontal separation in 
the perpendicular (cross-elevation) direction H. If the beam is pointed at the horizon at an azimuth 
of 0°, then V corresponds to elevation and H to azimuth. We define H - cos (p and V = -0 sin (p, 
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where 6 is the angular distance from the beam center and (f> is the azimuthal angle around the beam, 
with = 0° orthogonal to the antenna symmetry plane and = 90° corresponding to the downward 
direction in elevation (away from the GBT arm; see FiglUand Appendix lAll. 

The sunscans cover the most of the region -20° < V < +55° and -17° < H < +55°. The 
left-right reflection symmetry of the telescope implies that a sunscan at -H should be identical to 
one at H. The data were recorded at intervals of between 2 and 10 s, corresponding to separations 
of a few arcminutes on the sky at the adopted scan rates. 

The initial data set consisted of "vertical" raster scans spaced 2° - 4° apart. Among these, the 
"distant vertical" scans lying beyond H ~ 35° from the main beam (dotted lines at right in Fig.|6|l 
showed no evidence of sidelobes above the noise level. The theoretical values for these scans are 
small enough to be consistent with a measured value of zero, considering the uncertainties in the 
data. 

The second data set consisted of "horizontal" scans spaced by about 1° in V. The third set 
consisted of "vertical" scans, most spaced by about 1° in //, but including a group with 30' spacing 
refining the coverage of the range \H\ < 3°. Although we were not able to measure the power 
pattern to large negative angles H, both the theoretical calculations and our observations indicate 
that the pattern is symmetric about a vertical Une {H - 0°) through the center of the main beam. 

6.1 .2. Data reduction and relative calibration of the observations of the Sun 

A baseline was set using the longest scans, assuming that the sidelobes many tens of degrees from 
the main beam are relatively negligible; the calculations suggest that away from the spillover lobe 
the typical sidelobe has an amplitude of -75 dB with respect to the main beam. Backgrounds had 
to be subtracted, consisting of a large constant component plus a smaller elevation-dependent com- 
ponent with a slight amount of curvature which was well fitted by an exponential which accounts 
for atmospheric emission variations with elevation. The most distant vertical scans that showed 
no strong sidelobes were used to examine the elevation dependence of the background signal in 
the recorded data. Both background components varied to some extent from scan to scan, with the 
constant component also being slightly different for the two polarizations. Spikes in the data due to 
RFI were also removed. The numerous scan crossings were used to fix the amplitude of the shorter 
scans. Consistency of these scan crossings, and of the scans symmetric about H - 0°, indicated an 
accuracy of approximately 10% in the background- subtracted scans; the horizontal scans appeared 
somewhat less accurate, possibly due to the variation in the height of the horizon along a horizontal 
scan. 

At 1 .4 GHz the Sun can have significant temporal variations in its emission. Relative normal- 
izations of the three datasets were determined from the numerous scan crossing points. We could 
not find 21 cm measurements of the solar flux for the periods of our observations, and so we used 
the solar flux monitor measurements at 10.6 cm from the DRAc| to check that the relative normal- 
ization factors were consistent with the ratios of the solar flux. The second and third datasets were 
obtained near the solar activity minimum, and required only a small relative normalization factor; 
the first dataset had been obtained at at time when the solar emission was approximately twice as 
large, requiring a relative normahzation factor of about two. 

^ http://www.spaceweather.gc.ca/sx- 11-eng.php 
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Fig. 7. Estimating the optimal value of t/sI using three repeated observations at more than 4 x 
10'* positions near the NEP, each repeat having different stray radiation. For different assumed 
/7si, the difference between each observation of W and the mean of the three was calculated for 
every position. Histograms of the resulting distributions of differences show the sensitivity to 7751, 
including the dotted curve for no stray radiation correction. Correction is clearly beneficial, and 
a minimum dispersion occurs for ?7si ~ 0.097. Note that the dispersion does not quite reach the 
minimum predicted from line noise and baseline uncertainties alone (dashed line, see Sect. l7.2l i. 



Scans were smoothed along the scan direction, and linear interpolation was used to esti- 
mate values between the scans in a direction approximately perpendicular to the scan direction. 
Observations at negative azimuth angles were "flipped" to positive azimuth, providing additional 
coverage and consistency checks. Where scans crossed or approached closer than 30' to each other, 
a weighted average was used at their positions to avoid any sudden jumps. Separate interpolations 
were made for each of the three datasets, with a weighting approximately proportional to the den- 
sity of the observations on the sky. 

The result is a map of the relative GBT beam pattern covering the important sidelobes within 
the range \H\ < 32° (horizontal extent) and -20° <V < +55° (vertical extent). 



6.1 .3. Determining tine amplitude scale of the measured sidelobes 

The power pattern measured using the Sun needs to be scaled by an amount determined below, 
because the absolute brightness of the Sun at the times of observation is not known to the required 
precision. Because the sidelobes are fixed with respect to the GBT, as a given direction is observed 
the sidelobes wiU cross regions of different 21 -cm causing the stray radiation component to vary 
throughout the day. There are also different Doppler shifts. We took advantage of these temporal 
changes to estimate the scaling of the measured sidelobes by mapping a large area around the 
North Ecliptic Pole (NEP) containing > 4 x lO'* unique positions. HI spectra were measured at 
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each position on three separate occasions at different ranges in azimuth and elevation. For each 
position the spectra were reduced, calibrated, and corrected for stray radiation using Eq.|4]with a 
range of scalings which, when added to the calculated part of the sidelobe pattern, gave a range of 
values of rjsi- The amount of the telescope response in the calculated portion of the far sidelobes, 
which is generally the region alO > 60°, is 0.0187. 

A histogram of the differences between each of the three observations of the integral over the H I 
line profile (W, in K kms"^) and their mean is shown in Fig.Qfor different representative rjsi. The 
black dotted curve shows the large scatter in the measurements when there is no correction for stray 
radiation. The stray radiation correction is quite effective: the dispersion decreases to a minimum as 
the scaling increases, then subsequently increases. The optimum scaling of the empirical sidelobe 
pattern is the one that minimizes the dispersion, near t/sI = 0.097. 

A similar analysis, looking directly at differences in channel by channel instead of differ- 
ences in the area W, was perfor med on the NEP fiel d and on other areas that were observed multiple 



times during the course of the iMartin et alj ( 1201 ih surveys. For the spectra in the NEP field, this 



analysis yields a value of 7751 = 0.1014 ± 0.0018, where the uncertainty is solely a measure of the 
statistical accuracy of t^si; it does not take into account the accuracy of the sidelobe pattern. This 
places the above NEP W result within 2cr, indicating the consistency of the two analyses. For all 
additional (non-NEP) spectra - most of which only have two repeat visits - the T^b analysis yields 
a value of rj^i = 0.0947 + 0.0032. The statistical error was obtained by a chi-squared analysis; 
to determine the number of degrees of freedom, it was assumed that each continuously-observed 
"chunk" of sky (i.e., wherein all spectra have similar sky positions and LST values) could be con- 
sidered an independent measurement of t/si. There are proportionally a larger number of spectra 
in the NEP field, but they are all representative of a single patch of the sky with three repeated 
observations, and statistically will have more similar LST values — observations of nearly the 
same patch of sky at nearly the same LST do not yield entirely independent measurements of 7751. 
Obtaining the best %\ separately for all nine subregions observed in the NEP survey yields values 
that are all within Q.Scr of the best overall NEP value, according to their individual statistical errors, 
unlike the 13 separate non-NEP regions, where the scatter is consistent with the statistical errors. 
Thus the formal statistical error from the NEP measurement underestimates its true error, which 
is probably closer to the non-NEP error estimate, and so we decided that an average of the two 
values (NEP and non-NEP) would yield a less biased result for rj^j. Consequently, we have adopted 
77si - 0.0981 + 0.0023 to set the amplitude of the measured sidelobes on a correct scale relative to 
the GBT main beam. 



Note that ?/n,b + ?/si ~ 0.98; we cannot account for about 2% of the power. Cases were tested 
where an isotropic component rjiso was added to the sidelobes to bring the total beam power closer 
to unity. However, even rji^o - 0.01 (i.e., an increase of 0.005 in the sidelobe power seen by the 
sky) yielded significantly more overcorrection, with parts of stray-corrected spectra going nega- 
tive. The formal uncertainty of 0.0023 refers to overall increases or decreases in the measured 
sidelobe. Modifications to the shape of the measured sidelobes (e.g., an isotropic component; see 
also Appendix lB.6l l. would yield a different "best t/si" value. Considering possible variations yields 
a somewhat larger error estimate of 0.005. 
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6.2. Combining tiie calculated and measured power pattern 

The adopted power pattern P is a combination of the calculations of Sect. 15. II with the measure- 
ments of the Sun of Sect. 16. II Measured values atO > l°.2 begin to differ from the calculated values, 
while at < 0?8 the coarse, half-degree spacing of the measurements is insufficient to probe the 
smaller-scale beam structure there; for both these regions, differences between measurement and 
calculation can exceed a factor of 2. However, for 9 x 1°, measured values {P convolved with 
the Sun) typically agree to better than 20% with smoothed calculated values (P convolved with a 
paraboloid to yield a resolution of CK5) — this is nearly the best that can be expected from the inher- 
ent uncertainties in the measurements. Therefore, for CP8 < < 1°2, we make a smooth switchover 
from the (convolved) calculated P to the measured one. As we correct for sidelobes only at > 1°, 
this switchover has a negligible effect on the results. 

The theoretical values for P outside the region measured using the Sun are small enough to be 
consistent with a measured value of zero, considering the uncertainties in the data. Some of the 
back sidelobes are simply not accessible, being below the horizon; they are also not used in our 
evaluation of the stray radiation (Sect. 16.4b . Therefore at angles beyond where we were able to 
probe directly, the theoretical calculations of P were again adopted with a smooth switchover at the 
outer (f.5 edge of the measured region to prevent any possible discontinuities in the adopted P. 



6.3. Properties of the adopted response pattern 

Figure[8]is a contour map of the inner part of the derived GBT power pattern probed with the Sun; 
this is G in dBi. The obvious symmetry about a vertical line through the center of the main beam 
(H = 0°) is by construction, as discussed above in Sect. 16.1.11 Other displays of the GBT sidelobe 
pattern are given in Figs.|9]and[T0l these are on a linear scale and are for P - AkG. 

Away from the main beam, the power pattern is dominated by the forward spillover sidelobe, 
the arc of radiation from the secondary feed spilled past the subreflector (see Fig.|9]). In symmetric 
antennas such spillov er sidelobes are symmetri cal about the main beam, but with the offset sub- 



reflector of the GBT dNorrod & Srikanth 



19961) it is centered roughly on the cone axis defining 



the secondary, displaced by about V - 12° above the main beam in elevation (it retains left-right 
symmetry but is not circularly symmetric). This spillover lobe results from near-field (Fresnel) 
diffraction of the feed illumination from the sharp edge of the subreflector, which can be thought of 
as a disk occulting the sky. Our measurements of the diameters of the main spillover lobe and the 
fainter rings outside it are consistent with this, given the GBT geometry (Appendix iBl). The peak 
level along the ridge of this main spillover lobe is quite low, about +5 dB above isotropic and there- 
fore about 57 dB below the main beam, but because of the large area it contributes substantially to 
T/si and the stray radiation. 

The gap in the spillover lobe at {H, V) ~ (0°, 30°) is caused by the arm that supports the GBT 
subreflector, specifically a reflecting screen attached to that arm (Fig. 3). It deflects the spillover 
radiation back into the main dish over a 40° segment of the spillover lobe. The reflected radiation 
emerges on the sky as a pair of sidelobes well away from the main beam on the opposite side, 
at about H - +l°,y = -3° (Fig. [8]). Each peak is elongated by about 1°, approximately along 
the azimuthal direction as seen at the left in Fig. [10] They have a peak amplitude about ten times 
that of the spillover lobe, but comprise only a small part of the total beam integral, comparable to 
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H (deg) 

Fig. 8. Contours of the measured GBT far sidelobe gain G, relative to isotropic. The main beam is 
at the origin (indicated by the solid dot), and would peak at +61 .6 dBi. The H and V coordinates are 
the same as in Fig.|6] Elevation increases upward in the plot and the direction to the center of the 
subreflector is at V =12.3°. The main "spillover lobe" is outlined by the yellow dBi contour, with 
the blue 4 dBi and magenta 5 dBi contours defining its ridge; curving as a ring from below the main 
beam at (H, V) ~ (0°, -7°), it extends to about (H, V) ~ (12°, 28°), beyond which it is blocked by 
the screen on the telescope arm (note that this blockage also yields a ridge stretching upward from 
the subreflector center along the edge of the gap at an angle about 35° from the vertical; this ridge 
passes through the end of the spillover lobe). Complementary to this missing part of the spillover 
lobe caused by the screen are the twin peaks below the main beam, at (H, V) ~ (+1°, -3°), visible 
in the heavy black 10 dBi and thin gray 15 dBi contours. Outside the spillover lobe are three lower- 
amplitude rings (at radii ~ 26°, 29.5°, and 33° from the subreflector center). Inside is the Arago 
spot and surrounding rings centred at 11.879° on the axis above the main beam. Some features 
are easier to see in the 3-D plots of Figs. |9]and[T0l Outside the faint spillover rings and outside the 
region {\H\ < 35°, 20° < V < 53°), the sidelobe levels are taken from the antenna code calculations. 
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beam 




Fig. 9. GBT response pattern P, as measured from scans of the Sun, on a linear scale. The H and V 
coordinates are the same as in Figs.|6]and|8] Note that this plot has been rotated so that its features 
are more easily visible — the vertical coordinate V increases from left to right in this figure. The 
truncated peak at center-left in this figure comprises both the main beam at (H, V) = (0, 0) and the 
double-peak feature at (H, V) ~ (+1°, -3°); the peak of the main beam would he far offscale at 
P - 1.12 X 10''. Near the center of the spillover lobe, the Arago spot is visible. The gap at high V 
in the spillover lobe and its surrounding rings arises from the presence of the reflecting screen at 
the subreflector edge. 



the portion of the major ring removed in the wedge, qualitatively consistent with conservation of 
energy. This is discussed in more detail in Sect. IB. 41 

A more minor, but interesting, feature in the GBT beam pattern is the Poisson-Arago spot at 
H - Q°,V =11 .879° on the axis above the main beam (FigsiSl-fTOli. As expected from simple near- 
field diffraction theory, this is nearly the same direction as the center of the subreflector (12.3°). 
The width of the Arago spot and the set of rings seen around it are also in detailed agreement with 
simple diffraction theory (Appendix lB.2b . 

There is a feature in the (calculated) beam pattern located atV ~ -96°, with a peak amplitude 
about 57 dB below the main beam, slightly above isotropic (see Fig. lA.4l i. It has a width of several 
degrees in V, is somewhat wider in H, and is part of a ring, very asymmetric in both amplitude 
and position about the main beam. This sidelobe arises from spillover past the edges of the main 
telescope reflector (Appendix lAl). Along with this is another Arago spot in th e direction from the 
prime focus along the cone axis defining the primary dNorrod & Srikanthll996h . These "backlobes" 
contain roughly 2% of the total power, i.e., roughly 20% of the total sidelobes. However, because 
all but a tiny fraction of these are always below the horizon, they do not "see" the sky and so do 
not contribute significantly to the 21cm stray radiation. 
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Fig, 10. Measured GBT response pattern P, as in Fig.|9] but for sidelobes closer to the main beam. 
The Arago spot, peaking at s; 0.25, is visible at the right, along with its surrounding rings. The 
double-peak feature at (H, V) ~ (±1°, -3°) appears just to the left of the (truncated) main beam, 
with height P » 2.6. This contains the power scattered from the forward spillover lobe by the 
reflecting screen on the feed arm (Fig. |3]l. The peak of the main beam would be far offscale at 
P = 1.12 X 10\ 



6.4. Implementation of the stray radiation correction 

Stray radiation was calculated for the GBT spectra following the integral in Eq. [H with a program 
described more extensively in AppendixICl A model H I sky was constructed from the LAB survey 
data to give the input T\,(9, (p) for the integral, on a tiled grid in Galactic latitude and longitude. 
Given a specific GBT observation in a particular direction at a specific time, the GBT response was 
calculated for each tile location in the model sky more than 1° from the main beam. The model 
sky spectrum was accumulated from all directions above the local horizon after weighting by the 
beam response, accounting for atmospheric attenuation, and making the appropriate velocity shift. 
This integrated stray radiation spectrum was then subtracted from the observed T.^ spectrum, before 
correcting for atmospheric attenuation, and scaling by rj^b. Implemented in the language C on a 
modern workstation, the program can calculate the stray-radiation correction at a rate of several 
spectra per second, or of order a square degree per minute for our mapped data cubes. The program 
is now available for use by any GBT observer. 
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Fig. 11. Example of correction for stray radiation for a subregion in the NEP field. Upper panel: 
Dashed curves show spectra taken at different LST (black and gray) with different stray radiation. 
Each curve is the average of 280 contiguous HI spectra on the T^^b scale; averaging is essential 
for lowering the noise, to reveal more clearly the effects of stray radiation. Dotted curves: our 
calculation of the expected stray radiation. Solid curves: spectra corrected for stray radiation, now 
well aligned. Middle panel: mean difference of the corrected spectra (note that the difference of the 
uncorrected spectra, or of the stray radiation, would be far off scale). Lower panel: spectra showing 
the rms of the 280 individual differences before (dashed) and after (solid) the correction for the 
stray radiation. 

6.5. Examples of the effects of stray radiation 

Examples of the stray radiation and our ability to remove it are shown in Figs.[TT]and[T2l Depending 
largely on where the spillover lobe is positioned on the sky, the stray radiation correction ranges 
in peak amplitude from a fraction of a K to several K. Most stray radiation tends to lie relatively 
close to zero velocity but is typically broader in velocity than the corrected/intrinsic spectra; the 
total velocity range over which stray radiation exceeds 0.1 K varies from a few tens of kms"' to 
about 150 km s"'. 

In these examples we compare spectra of the same region taken at different LST so that they 
will have different stray radiation. To reduce the noise and reveal the subtle changes after correction 
for stray radiation we have averaged many contiguous spectra that will have rather similar stray 
radiation. Figure[TT]illustrates a case where there are large stray corrections that produce consistent 
results. For one of the epochs (gray lines in Fig.fTTI). the stray correction removes excess emission 
on both wings of the main peak, and the stray correction is significant all the way out to v ~ 
-120 kms Compared to the up to 2 K differences in the uncorrected spectra, the difference 
between the corrected spectra is an order of magnitude smaller (middle panel). The difference after 
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Fig. 12. Like Fig. [TT] but based on 264 contiguous spectra in a subregion of Nl. The stray cor- 
rection at one LST is slightly too large, pushing the corrected spectrum slightly negative in the 
high-velocity wing. However, overall the correction is as effective as in the NEP example; note that 
the middle and lower panels are on the same scales as in Fig.fTTI 

correction is not, however, zero, and some residual effects of stray radiation must still be present in 
the spectra. 

Figure [12] provides another illustration from the Nl region, a region of low signal where the 
stray correction can have a relatively large effect on the derived line profile and its integral W. Here, 
for one observation (gray Unes in Fig.[T2ll the stray correction overcorrects somewhat on either side 

of the main peak, yielding a slight dip below zero near v ~ 30 km s ' and v 60 km s ' . However, 

the stray correction still significantly reduces the difference between the two observations and the 
amplitude of the rms difference spectrum, as shown in the middle and lower panels of Fig. fT2l 
respectively. Note that these panels are on the same scales as the corresponding panels in Fig.fTTI 
the residual effects remaining in the corrected spectra are similar. 

The quantification of the errors remaining because of imperfections in the stray radiation cor- 
rection, and the cumulative effects on W, will be discussed further in Sect. l7.3l in the context of the 
other sources of error. Judging from many more comparisons of repeat measurements of different 
fields, the largest errors in the stray radiation correction tend to occur within a few tens of kms ' 
of zero velocity, with a typical amplitude of 0.1 to 0.2 K, although in the worst cases the error can 
exceed 0.5 K. 

7. Error estimates 

In this section we evaluate the various contributions to the errors in the GBT H I spectra. The effects 
of errors are manifested in the non-reproducibility of measurements. Comparisons between spectra 
can be done on a channel by channel basis, T{v), or for a line integral W (in K km s ') of T(v) over 
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some velocity range. Because of different applications of HI spectra, it is relevant to address the 
errors for each metric. 



7. 1 . Line noise 



The individual survey spectra from lMartin et al 



( 201 1 ) have an rms noise in emission-free channels 



of ctq = 0.16 K. When th e spectra are interpol ated into the data cube at Nyquist sampling this is 



reduced to cry - 0.11 K (iMangum et al.l 120071) . The 21-cm line emission itself can significantly 



increase the noise at velocities where it is bright: 

cr(v) = cro(l + T(v)/T,y,), (6) 

where T^y^ is the system temperature (Tsys ~ 20 K, Sect.|2|. For clarity, T(v) is referred to simply 
as T for the remainder of this discussion. 

All the GBT H I observations measure spectra in two orthogonal linear polarizations, labeled 
XX and YY, which should be receiving essentially identical 21-cm emission and identical stray 
radiation. Therefore, we are able to check the above equation by comparing spectra in the two 
polarizations. Both were processed in identical parallel streams prior to removal of distinct instru- 
mental baselines. The difference in the baseline-subtracted XX and YY spectra provides a good 
indication of the error inherent in a single measurement. It eliminates the additional uncertainty 
arising from the stray radiation subtraction which is in common, but includes channel noise and 
baseline error. This difference would also reveal any mis-calibration of the two receiver channels 
and any real differences in the received H I signal because of differences in beam shape, pointing, or 
sidelobes. These latter effects, however, are thought to be small compared to the other error terms. 

Because the errors in XX and YY are independent, and ultimately the XX and YY spectra 
are averaged together to form T = {Txx + Tyy)/!, the estimator of interest in assessing errors 
in T is the dispersion o-j, the standard deviation about the mean of AT^poi - {Txx - Tyy)/2. 
Figure [T3] displays this dispersion based on all of our NEP spectra. The data over the v range 
of -50 kms^'to +25 km s 'have been binned by T, each bin containing 2 x 10^ points. The T- 
dependence of the standard deviation follows the prediction from Eq. |6l overplotted in Fig. [13] 
using (To = -^(o"g^)^ -I- {(Tq^)^/2 = 0.111 K as measured in the individual spectra and a typical 
T,y, of 20 K. 



7.2. Baselines 

The mean of the distribution of AFpoi within each bin, also shown in Fig. [T3] is slightly offset 
from zero, indicating a systematic difference in XX and YY that can be attributed to imperfect 
baseline removal. This offset - though well-characterized by a mean of -15 mK - does vary from 
spectrum to spectrum with a standard deviation of ~ 27 mK. This gives some sense of the size of 
the baseline error in the polarization-averaged spectrum T; it is so small compared to cr(v) that it 
has little effect on the dispersion described above. However, because a baseline error is systematic 
over many channels, it can accumulate as a significant error in W. 

We illustrate this using the same NEP data, defining W to be the line integral over A^ch - 94 
channels of width Av = 0.8 km s ' over the range -50 km s ' to H-25 km s ' . We calculated the dis- 
persion (Tw of AWpol = (Wxx - Wyy)/2 for bins in W (2400 points each), plotting this in the upper 
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Fig. 13. Estimating the noise in a spectrum T using ATpd. Open circles give the standard deviation 
of AT^poi within bins of 2 x 10^ data points as a function of T, using spectra from all NEP pointings. 
This agrees closely with the prediction (dotted line) from Eq.|6]with ctq = 0. Ill K and T^y^ =20 K. 
Also shown is the mean of AFpoi for each bin. 



part of Fig. [141 The upper dashed line is a prediction of this dispersion, combining in quadrature 
the minimal channel noise 0.111 V-^chAv (dotted line) and a baseline error of 0.027A^(;hAv (both in 
K km s '); note that the latter is now the larger contribution, because of how it accumulates system- 
atically. The larger values of W often result from larger T and thus larger cr(v), although this is not 
necessarily the case with very broad lines. The slight upward trend in the observed dispersion to- 
ward larger W could therefore be a result of the increasing contribution of o-(v) to the overall error. 
For completeness, we show the mean of AWpoi as well and the prediction -O.OlSA^chAv K km s ' . 

The uncertainty in W = (Wxx + VVyy)/2 should also be on the order of 2 K kms for opti- 
cally thin emission, this corresponds to an uncertainty of only 4 x 10'^ cm"^ in column density. 
Note, however, how this varies with the number of channels used in the W integral; it should be 
determined self-consistently for the W appropriate to each different region or application. 

An independent estimate of the baseline errors can be made directly from the third-order poly- 
nomial models used to fit the residual baselines in the spectra (13.31 ). A Monte Carlo analysis 
based on the uncertainties of the coefficients of the Legendre polynomials yields errors in W of 
~ 0.7 K kms"' over the same -50 kms"' to -i-25 kms"' velocity range, suggesting that the errors 
derived from A Wpoi are an upper limit. For the remainder of this paper we adopt this upper limit as 
it includes potential systematic errors (e.g., offsets) between XX and YY which are not detectable 
in their average, W. 

The good agreement between the data and the predictions in Figs. [T3] and [14] indicates that 
we have a good understanding of the origins of the errors that arise from noise and instrumental 
baselines and that there are not large differences in the H I signal measured in the two polarizations 
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Fig. 14. Estimating the noise in the Hne integral W using AWpoi- Open circles give the standard 
deviation of AWpoi within bins of 2400 data points as function of W, using spectra from all NEP 
pointings as in Fig. [T3] This agrees closely with the prediction (dashed line) based on the accu- 
mulation of line noise (dotted line) and baseline errors. Note that the baseline error is dominant in 
measurements of W, even though it is not dominant for measurements of the spectrum T (Fig.fTSl). 
Also shown are the data (filled circles) and the prediction (lower dashed line) for the mean AWpoi. 

of a single observation. However, we still need to assess the errors arising from the stray radiation 
correction (Sect. 1731 ). and any other time-varying error contribution (Sect. lT4] l. 

7.3. Stray radiation 

As seen in the rms curves in Figs.[TT]and[T2] the baseline error and line noise are not the entire story. 
Changes that affect both XX and YY simultaneously and systematically can only be diagnosed with 
repeated observations at the same position. 

The data from the NEP field, covering 44100 spatial pixels three times each, are used here to 
examine the reproducibility of spectra with a large range in T and W. This allows us to assess 
how the uncertainty in the stray radiation correction contributes to the overall error. To estimate 
the errors in a single observation, T = {Txx + 7^yy)/2, for comparison with the results in the 
subsections above, we examine the statistics of ATij - {Ti - Tj)/ V2, where / and j denote two 
separate observations, obtained at different LST and elevation as this region was mapped over 
several months (2006/10 to 2007/01). 

The dispersion of ATij for data binned (2 x 10^ per bin) in T, as in Sect. 17. II is plotted as open 
circles in the upper part of Fig.[T5]for each of the three observation pairs. The expected standard 
deviation from line noise (plus a minimal baseline component of 0.027 K added in quadrature) is 
plotted as a dashed line for Tsys = 20 K. Also shown as solid circles is the mean of ATij in each 
bin. We attribute both the excess of the observed dispersion above this prediction and the non-zero 
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Fig. 15. Estimating the eiTor in a spectrum T following stray radiation removal, using repeated 
observations in the NEP field. Open circles give the standard deviation of ATij within bins of 
2x10^ data points as a function of T, using spectra from the three pairs of NEP observations (dark 
gray: 3-2, light gray: 1-3, black: 2 - 1). These lie slightly in excess of the prediction based 
on only line noise plus baseline errors (dashed line). Also shown is the mean of AT/j for each bin 
for each ij combination. Both the excess and non-zero means arise from imperfect stray radiation 
removal. Crosses show a simple estimate of the total error which combines 6% of the rms spectrum 
of the stray radiation correction in quadrature with the line noise. 



means to errors or imperfections in the stray radiation correction. A simplistic estimate of the total 
error can be obtained (crosses in Fig. [TSl l by including some fraction (here 6%) of the rms of the 
stray radiation corrections, in quadrature with the line noise. 

The mean difference does vary from spectrum to spectrum with a standard deviation of 
~ 0.041 K. This is about 40% of the typical dispersion in T attributed to stray radiation. This 
gives some sense of the size of the systematic error from the stray radiation correction (see also 
Figs.fTTIandfTSb: it is small compared to the dispersion of ATij but because it can be systematic over 
many channels, it can accumulate as a significant error in W. This is illustrated using the same NEP 
data. Again defining W to be the line integral over Nch - 94 channels of width Av - 0.8 kms"^ 
over the range -50 km s"' to -h25 km s"' , we calculated AWij - (W,- - Wj)/ V2 for bins (800 points) 
in W. As shown in Fig. [16] the mean of AWij is typically < 1.5 K kms"' (lower points) across all 
W bins. The dispersion in AWij over all ij combinations is plotted as open circles in the upper part 
of Fig. Hi 

This analysis of observations of W in NEP, taken over a large range of time and with very differ- 
ent stray radiation corrections, indicate that the data are reproducible to an accuracy of 3 K km s ' 
in a single mapping. This represents an accuracy of 1% to 3%, depending on W. As in Fig. [T4l 
the dotted line in Fig.[T6]is the predicted dispersion from channel noise alone, whereas the dashed 
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Fig. 16. Estimating the noise in the Hne integral W using AW,j. Open circles give the standard 
deviation of AW,j over all ij combinations within bins of 800 data points as a function of W, using 
repeated observations of spectra from all NEP pointings as in Fig. [15] This agrees closely with the 
prediction (dashed-dotted line) based on the accumulation of systematic errors in the stray radiation 
and baseline corrections added in quadrature with the line noise. The line noise and baseline errors 
(dotted and dashed lines from Fig. HM are shown for comparison. Note that the error arising from 
the stray radiation and baseline correction is dominant in measurements of W, even though it is not 
dominant for measurements of the spectrum T (Fig.fTSl). Also shown are the data (filled circles) for 
the means AWjj for each ij combination. 



line accounts for baseline errors as well. The upper dashed-dotted line is our prediction combin- 
ing the channel noise in quadrature with a systematic error from the stray radiation correction of 
0.041A^chAv. This systematic error includes both baseline and stray radiation errors as it is not pos- 
sible to disentangle these here. The prediction is in close agreement with the observed dispersion. 

The largest source of uncertainty in the stray radiation correction appears to be from incom- 
plete knowledge of the sidelobe pattern for the telescope. Other sources of uncertainty in the stray 
radiation correction are discussed in Appendix iDl Note that the velocity range used in this exam- 
ple was selected to accentuate errors from stray radiation. Because there is little significant stray 
contamination at \vlsr\ » kms measurements of W for intermediate and high-velocity HI 
components will have uncertainties closer to those predicted from the appropriate AWpoi. 

For high-latitude regions of very low and W, the stray radiation from sidelobes overlapping 
H I near the Galactic plane can be as strong as the actual signal from the main beam at some 
velocities. The highest accuracy measurements can be obtained by timing the observations so that 
the spillover lobe does not lie near the bright H I in the Galactic plane, since that minimizes the 
stray radiation and the errors associated with its removal. With our knowledge of the extent and 
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location of the sidelobes, we are able to d etermine the ideal L ST range and successfully apply this 
strategy to the very faint ELAIS Nl field ( Martin etallboil ). 



7.4. Repeated observations of H I calibration standards 

The long-term reproducibili ty of the GBT s pectra can be gauged usi ng our observations o f standard 



HI caUbration regions, S6 (lWilliamJl973L hereafter W73) and S8 jKalberla et al.ll 19821 hereafter 
KMR), made over the period 2005/10 to 2008/03. The integration time for each spectrum obtained 
was 180 s and so these have much lower Tsys noise than the typical map spectra. Nevertheless they 
are affected by other sources of error, principally in the stray radiation correction. The spectra were 
reduced using the standard procedures (Sect. |4|. For each region we created an average spectrum 
(T) from all the data. 

There were 31 and 29 repeated observations of S6 and S8, respectively, each consisting of two 
spectra, Ti and T2, from consecutive 180 s integrations. The stray radiation correction, and also 
the baselines, ought to be very similar for the two spectra in each pair. We therefore examined 
AT12 - (Ti - To)/ V2 channel by channel. The standard deviation of AT12 for the 29 pairs in S8 
is shown by the plus symbols in the lower part of Fig.[T7] for channels in the velocity range used 
for the KMR calibration (Sect. 18. 11 1. As expected, this tracks closely the prediction (dashed line) 
based on line noise alone from Eq.|6]with ctq - 0.026 K as measured in emission-free channels, 
and T,y, =20 K. 

As a second check, we examined the standard deviation of AT^poi as in Fig.[T3](Sect. l7.1l l. These 
data are drawn as open circles in Fig.[T7]and lie slightly above the dotted line because of baseline 
errors, which for these long integrations with lower Tsys noise have relatively more importance. 

Next we examined the standard deviation of AT - T - {T}. This is shown by the upper filled 
circles in Fig. [17] This somewhat larger dispersion (but note that it is still less than 1%) is from 
additional errors from the stray radiation corrections and potentially tiny changes of gain with time. 
We computed the rms spectrum of the actual stray radiation corrections and take some percentage 
of this as a rough estimate of what the error of the stray radiation correction might be. The crosses 
in Fig. [17] show the result of taking just 7% of the rms spectrum and adding this error estimate 
in quadrature to the other errors given by the standard deviation of ATpoi- Although a simplistic 
description, it provides a reasonable explanation for both the magnitude of the standard deviation 
of AT and its dependence on (T). 

Pursuing another approach that might reveal, for example, gain variations, we looked at the 
channel by channel correlation of T with {T}, the slope of the regression being s. We carried out 
these regressions for data in the velocity ranges as defined by W73 and KMR for their calibration. 
By definition s will have a mean of unity, but it might vary from observation to observation due 
to various errors. Histograms of s are presented in Fig. [18] The standard deviations are 0.008 and 
0.006 for S6 and S8, respectively. The GBT spectra are thus reproducible to better than 1% and thus 
any changes in gain are smaller than 1 % as well. The smallest values of s occurred for observations 
taken during periods of rain, suggesting an additional source of atmospheric opacity not currently 
included in our reduction procedure; these outliers (one pair for S6 and two pairs for S8) were 
omitted in the above analysis of AT". 
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Fig. 17. Estimating the noise in the spectrum T of S8, using repeated observations. Plus symbols 
show the observed standard deviation of AT12 from channel by channel comparisons of two con- 
secutive spectra, which ought to have very similar stray radiation and basehnes. The measurements 
agree closely with the prediction (dotted line) from Eq.|6]for line noise alone. Open circles give the 
standard deviation of Arpoi; these lie slightly above the prediction because of baseline errors. The 
filled circles give the standard deviation of AT - T - (T), which reflects all errors. As a prediction 
of this, the crosses show the result of adding 7% of the rms spectrum of the stray radiation cor- 
rection in quadrature to the standard deviation of Arpoi. The diamonds result from further addition 
of the effects of tiny 0.3% gain changes. The S8 data are consistent with errors in stray ~ 7% and 
scaHng (e.g., gain) errors of < 0.3%. 



These two analyses of S8 are not decisive as to the presence of tiny secular changes in gain 
along with the errors from the stray radiation correction. S6 provides an opportunity to discriminate 
because the stray radiation is relatively weaker compared to the T^b signal. Repeating the same 
analysis as for Fig.[T7]for the case of S6, we found that the dispersion in AT", though smaller than 
in Fig. [17] is too large and has the wrong {T} dependence to be explained by an error of 5 - 10% 
of the rms stray radiation correction, suggesting that gain changes of order 0.5% might be present. 
In summary, a common explanation for S6 and SB suggests gain changes of order 0.3% to 0.5% 
combined with stray correction error of size 7% of the rms stray radiation correction. The additional 
effect of 0.3% gain changes in S8, added in quadrature with the other errors, is illustrated by the 
diamonds in Fig. [T7] We have also reproduced Fig. [T8]for the case of S8 using simulated spectra 
starting with {T} and adding line noise, 7% of the rms stray radiation correction, and 0.3% gain 
changes with independent random seeds, verifying the self-consistency of this identification and 
quantification of the sources of error. Note that in the NEP (analysis in Fig.fTSll and our other HI 
survey regions the signal T^b is so small that any tiny gain errors would have a small effect. 
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Fig. 18. The distribution of slopes s from regressions of T on (T) for S6 (dashed) and S8 (soHd) 
from observations made over 2.5 years. The dispersions of 0.007 and 0.008, respectively, result 
from line noise, baseline errors, and errors in the stray radiation correction, and potentially from 
gain variations. 



Robishaw & HeilesI (120091) 



This excellent reproducibility is to be contrasted to the assertions by 
that the GBT spectra suffer from 10% errors in the calibrated gain. Their analysis depended on 
solving simultaneously for gain and several parameters describing the amplitude of their adopted 
sidelobe pattern. However, their sidelobe pattern was only approximate, so that the resulting rela- 
tively large errors in the stray radiation correction become confused as implying significant changes 
in gain. We have seen such an ambiguity above, but at a level < 1%, and our joint analysis for S6 
and SB via Fig. [17] and Fig. [18] is consistent with contributions from errors in the stray radiation 
correction and only tiny secular gain changes. The tiny gain changes found are consistent with what 
might be expected from the effects of gravitational and thermal distortions on the passive surface 
of the main paraboloid of the GBT, as described in Sect. [2] 



7.5. Error budget and benefits from repeated observations 

The contributions to the total error in the HI spectral data are summarized in Table [3] for both the 
T spectra and the integrated W. In the case that these errors are all independent, these can be added 
in quadrature to assess the total error. 

Observations are often repeated to "beat down the errors". This is certainly beneficial for the 
errors arising from line noise and even for the baseline errors. However, it is apparent that the dom- 
inant source of error in W for our NEP map is from imperfect stray radiation removal. Furthermore, 
it is not clear how the dispersion of AW,/ is to be interpreted in assessing the errors in even a single 
observation. In the worst case, the measured dispersion could be entirely due to the uncertainty 
in one of the observations, and so the uncertainty in that single observation of W would be V2 
times larger, on the order of 4 K kms ' or up to a few percent of the values typical in NEP. An 
error of this size corresponds to roughly 10% of the typical stray radiation contamination for these 
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Table 3. Summary of enors in GBT H I spectral data 



Error type 


cr" of T 


Error in W'' , accumulated 




(K) 


from o" of r (K km s"' ) 


Tsys noise 


0.25'- 




line noise, cr(v) 


/3V(4 sVrd + rM/r,,,)" 


VAU rms(cr(v))Ai' 


baseline 


0.027 


0.027AfchAv 


stray 


o.07r,t,,,(v) 


0.07W,„ay 


scale*" 


o.oo5r(v) 


0.005W 



Notes. Except where otherwise noted, all cr are on the T^f, scale, adjusted for 77,^1, = 0.88. Error 
calculated for A'ch channels with spacing Av. Assuming that each W error is independent, the total W error 
is the sum of these, in quadrature. RMS noise in antenna temperature (T.J for 1 second of integration in a 
1 km s"' channel for two polarizations and in-band frequency switching. RMS noise of Tmb in a 0.8 km s"' 
channel for two polarizations. /? = 0.16 for a single pointing and /? = 0.11 for spectrum interpolated into the 
data cube. Integration times are t = 4s for our scanning strategy (see Sect.|3]l. For long S6 and S8 integrations, 
t = 180 s (see Sect. l7.4t . Tsys is approximately 20 K. Includes any gain variations. 

NEP observations, compared to the estimated 7% in Table[3] Note that additional observations will 
not necessarily make the resulting average more accurate. One can hope that there will be some 
cancellation, but it will depend on the accuracy of the stray radiation correction for the times of 
observation. 



8. Absolute calibration 

8.1. Standard H I calibration directions 

S8 at (l,b) - (207?00, -15?00) is among the lAU primary HI calibrators and has been studied 
extensively by KMR using the Effelsberg 100-m telescope. To compare our repeated GBT obser- 
vations with KMR we computed the line integral W over the prescribed velocity range -5.1 < v < 
+22.3 kms finding Wgbt = 831 ± 5 Kkms ', the standard deviation being consistent with the 
expectation from Fig. [T7] At the same angular resolution, KMR found VVEffei = 846 + 14 K km s"' , 
the estimated eiTor combining a general scale uncertainty of 1.5% and an approximate 0.2 K sys- 
tematic uncertainty from their stray radiation correction. Thus WEftei/WcBT - 1-018 + 0.018. We 
conclude that our independent calibration procedure produces 2 1 -cm spectra in agreement with this 
detailed previous calibration at the level of the uncertainties, 2%. 

S8 is a calibrator (along with the secondary calibrator S7) for the Leiden/Dwingeloo Survey 
which forms the northern p art of the LAB survey . The LAB resolution after regridding the spatially- 



sampled data is about 40' (IKalberla et al.ll2005b . To study changes in the spectrum with angular 



resolution, KMR made a small 1° map sampled on a rectangular grid at 5' intervals centred on 
S8. Similarly, we made a 1°.5 map with our standard scanning setup. Errors from the stray radiation 
correction and baselines ought to be fairly uniform over the map. Convolving our map to ever lower 
resolution, we reproduced the resolution dependence found by KMR (their Fig. 5). From the GBT 
resolution to the LAB resolution, W increases by a factor 1.017 + 0.009. Using this to scale our 
pointed observations gives W = 845+9 K km s"\ to be compared to 849+9Kkms ' independently 
from our convolved small map, and 856 ± 14 K km s ' found by KMR. The value computed from 
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the LAB cube (with the Wakker correction - see the following section) is 844 ± 12 K kms"', in 
close agreement at less than the 1% level. 

The other region that we measured repeatedly was S6 at (I, b) - (1?91,41°42). Determining W 
over the range -6.86 < v < +5.86 kms ' (approximating W73; see KMR for a discussion) gives 
Wgbt - 289 + 2 K kms"'. From a r.5 scanned map we find that the correction to the 35' (40') 
angular resolution of Hat Creek (LAB) is 0.990+0.006 (0.993+0.006) giving 286+3 K km s ' (287+ 
3). The value directly from the convolved map is 289 + 3 K kms '(290 + 3). S6 is not a primary 
standard and has not been measured as accurately as S8. W73 reports 299 + 22 K kms 'which is 
probably systematically several percent high because it has not been corrected for stray radiation. 
Within the large 7% uncertainty, it agrees with our measurement. The value computed from the 
LAB data, interpolated to the S6 position, is 292 + 5 K kms ', again in close agreement with our 
measurements at the 2% level. 



8.2. Comparison with the LAB survey 

For comparison with the LAB Survey, the GBT H I maps were convolved to the 40' angular res- 
olution of that survey. The convolution results in spectra with a negligible error from line noise, 
but any errors from the baseline and stray radiation corrections are not significantly reduced. For 
repeated observations of a region, we used the average spectrum. 

Where the GBT observations completely cover a LAB beam, we interpolated the convolved 
GBT data to the LAB positions in Galactic coordinates ((f.5 x 0°5 grid in I and b) and to the slightly 
coarser velocity grid. The LAB and convolved-GBT spectra almost always exhibited very similar 
spectral features, usually with a small difference in overall scale. The line integral W is used to 
quantify this scale; W is defined as the integral over all channels with T > 1 K in either spectrum. 

We performed a regression of Wlab on VVgbt, and since the LAB-survey-cube points at high 
Galactic latitudes b are not independent, weights of cos b were applied. This regression yields 

Wlab =(1.0298 + 0.0023)x Wgbt + (1-96 + 0.21) K kms-'. (7) 



Wakker et al.l (1201 In has recommended that a Gaussian of peak amplitude 0.048 K, FWHM - 
167 kms" '(equivalent to 8.5 K kms ' over the entire LAB velocity range), and center v - 
-22 kms ' be subtracted from all LAB spectra. After making this "Wakker correction" to the 
LAB spectra we find 

Wlab = (1.0248 + 0.0022) x Wgbt + (0.71 + 0.20)Kkms-'. (8) 

The intercept is now significantly lower, and so our data support the application of the Wakker 
correction to the LAB Survey data. The following analysis adopts this correction. However, note 
that this correction is statistically based and might be different in different regions of the survey 
and from spectrum to spectrum. 

If we adopt the hypothesis that the true intercept is zero, then 

Wlab = (1.0288 ± 0.0012) x Wgbt- (9) 

We conclude that overall the LAB scale is about 3% higher than for the GBT. This is taken into 
account in the model H I sky used for stray radiation correction; see Appendix lC.il 
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Fig. 19. Wlab/VVgbt for every spectrum that could be compared (dots, color-coded by region). 
Filled circles show, for each of the observed regions, the weighted average ratio and its uncertainty 
(smaller than the filled circle for the typical uncertainty of ~ 0.004 in the regional averages), along 
with the region's average Wgbt and its dispersion (indicated by the horizontal errorbars). These 
comparisons suggest that a single scale factor (horizontal solid line, from Eq. |9ll might be too 
simplistic. 



A single scale factor or a single intercept might be too simplistic. This is illustrated in Fig. [191 
which s hows the ratio VVla b/W^gbt for every spectrum compared. In agreement with the conclu- 
sions of iHiggs et al.l ( l2005b . we find that the LAB Survey appears to have random "calibration" 



errors, typically of a few percent. About 20% of the LAB spectra appear to be mis-scaled (up and 
down) by between 5% and 10% and another 7% of the spectra by more than 10%. There are a few 
very large outliers where a LAB spectrum has been mis-scaled by as much as 30%; these outliers 
were not included in the above analysis (though they are plotted in Fig. [T9l ). Errors in the LAB 
Survey are discussed in more detail in Appendix ID. 21 (note that random LAB calibration errors 
tend to cancel out in the stray radiation calculation). 

Also shown in Fig.[T9]are the weighted average ratios and their uncertainties for the 17 regions 
that we mapped, as obtained from regressions such as in Eq. |9] for each region. The diff^erences 
from region to region are significantly more than the formal errors, suggesting that there might be 
systematic normalization errors of a few percent in the LAB survey data which vary with position 
in the sky. However, we have not found any systematic trend in the Wlab/W'gbt ratio as a function 
of position. 

Fig. [19] does not show any clear trend of Wlab/W'gbt or its dispersion vs. W. The S6 and S8 
maps discussed above (Sect. 18.11 ) have the largest W values; based on six independent comparison 
positions in each of them, we find Wlab/Wgbt = 1.015 + 0.010 and 1.019 ± 0.007 for S6 and 
SB, respectively. Restricting the comparison to the four positions farthest from the GBT map edges 
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reduces the ratios to 1 .007 ± 0.008 and 1 .010 + 0.007, respectively. Though smaller than the overall 
ratio, these are still compatible with a 3% difference between LAB and GBT scales. 



8.3. Comparison with Lyman-a measurements of H I 



Using a sample of 28 bright QSOs and AGN as targets. iWakker et alj (1201 II) have compared mea- 
surements of the foreground Galactic A^hi as determined by UV spectroscopy in the Lyman-ff ab- 
sorption line with measurements of A^hi from 21 -cm GBT observations in the same directions. The 
GBT data were reduced as described in this paper, and the correction for stray radiation was often 
significant. Third-order polynomials were fit to channels between -300 < visr < - 1 50 km s ' and 
H-lOO < vlsr < +200 kms"'. For some directions showing emission from high-velocity clouds 
the lower velocity range was changed to -300 < vlsr < -200 kms"'. Each 21-cm spectrum was 
examined for quality and a third-order polynomial baseline was found to be a good fit to the instru- 
mental baseline. The ratio AHi(LyQ')/AHi(21 cm) - 1 .00 ± 0.07 (Icr) indicating excellent agreement 
between the two entirely independent sets of measurements and independent tracers of H I. Some 
of the scatter must certainly result from structure in H I within the GBT beam as there is an enor- 
mous difference between the angular scales sampled by the UV absorption and 21-cm emission. 
We take this result as confirmation of the accuracy of our overall calibration procedure and an 
indication that the total error in GBT values of Ahi must be significantly less than 7%. 



9. Summary and conclusions 



This paper describes the results of a program to develop a calibration procedure that allows ac- 
curate measurement of 21-cm HI spectra with the Green Bank Telescope of the NRAO. Using a 
combination of measurement and calculation we have developed a model for the all-sky response 
of the telescope and use it to correct for stray radiation in the GBT 21-cm spectra. 

Several methods were used to estimate errors in the final spectra. Stochastic noise and instru- 
mental baseline uncertainties are well understood and can be made quite low in long integrations 
(Sect. 17.1b . The overall calibration to an absolute brightness scale appears to be correct to within 
a few percent as judged against measurements of standard radio continuum sources and the moon 
(Sect. IS. 

The correction for stray radiation is the most uncertain, and can produce errors at any given 
velocity of as much as 0.5 K, though these occur mainly at IvisrI < 20 kms"' (Sect. l63] i. Errors in 
the total W caused by errors in the stray radiation correction are typically less than 3 K km s ' or 
4 K kms ' in the worst case (Sects. [73] and l73] l. These are equivalent to an error in optically thin 
Ahi of ~ 5 X 10'** and 7 x 10"* cm"^. Overall the system i s quite stable. We see no evidence for the 



'10% gain fluctuations" reported by 



Robishaw & Heiles 



( 2009 ). A n independent investigation of 



2007!) achieved accuracies of 3.5*5 



the HI content of nearby galaxies with the GBT (iHogg et al. 
where the main source of error was baseline uncertainties. 

The corrected GBT data are in good agreement with other measurements of H I. Most inter- 
estingly, the GBT data give, on average, identical va lues of A^m as those derived from Lyman-o- 

and our Sect. 18.31 1. The 



measurements towards a sample of AGN and QSOs (IWakker et al 



2011 



GBT data also agree to within a few percent with previous measurements of "standard" H 1 cal- 
ibration directions and with other H I observations at lower angular resolution. We have shown. 
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however, that applying the "Wakker correction" to the LAB survey JwakkeretalJl201lh improves 
the agreement between that survey and the GBT spectra. 

For the GBT at 1.42 GHz we find that 0.88 of the telescope's response is within 1° of the main 
beam, with most concentrated within (f.2 of the main beam, while another 0.098 ± 0.005 of the 
telescope's response is in sidelobes more than one degree from the main beam. We thus account for 
0.978 ± 0.005 of the telescope's response including that which always lies on the ground, although 
the uncertainty of ±0.005 applies only to that part of the sidelobe that sees the sky: the part that 
sees only the ground cannot be measured, only calculated. These results are very constrained. Given 
the measurements of the sidelobes there are only two parameters that can be varied: z/mb and t/sI- 
The value of //si has been optimized from the observations (Sect. l6T3] i: any change in Tj^i, would 
increase discrepancies with other measurements (Sect.|7]). Attempts to add an isotropic component 
to the GBT beam pattern gave unacceptable results for the stray correction, implying //iso < 0.01. 
Strictly speaking, only the part of the isotropic component lying above the horizon is constrained 
by this analysis. However, the magnitude of any extra component of P in a backlobe is also strongly 
constrained by the low receiver T^ys- Given that our model already has nearly 2% of P intersecting 
the ground, any additional component that lies on the ground must have much less than 1% of the 
total P. 

The cahbration techniques are illustrated using data from HI maps made in connec 



2011 



tion with studies of interstellar dust and the cosmic infrared background (IMartin et al. 
Planck Collaboration et alj[201 la ). Overall, the data reduction process described here reduces sys 



tematic uncertainties in GBT H I spectra by at least an order of magnitude. Extra precision might 
be obtained with a better model for the far sidelobes including the possibility of reflections from 
the ground, allowance for the slight changes in telescope geometry with elevation angle, and more 
detailed consideration of atmospheric opacity during periods of rain. Given the unblocked optics 
of the GBT it should be possible, in principle, to construct a receiver for 21 -cm work with negli- 
gible forward sidelobes and a stray radiation component tha t is likewise negligible. This could b e 



achieved in the near future with phased-array feed receivers (ILandon et alJl201C ; 



■Teffset al.ll2010h . 
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Fig. A.l. Geometry of the GBT used in Antenna Workbench calculations. SF and PF are the two 
foci of the ellipsoid of which the subreflector is a segment. PF is also the primary focus of the 
main parabolic reflector. The "s" coordinate system used for the near-field subreflector patterns is 
aligned with the Zs axis running along the ellipsoid axis from PF toward SF. The illumination of the 
primary by these patterns is given in the "/?" coordinate system, rotated with respect to "i", so that 
the Zp axis runs from PF toward the projected center of the primary. The all-sky far-field response 
is given in the non-subscripted coordinate system, where the z axis runs from PF in the direction of 
the main beam. All have y and z in the symmetric plane of the antenna. 



Appendix A: Calculation of the GBT antenna pattern 



The main beam and the sidelobes of the GBT at 1.40 GHz were calculated using a reflector 
antenna code developed at the Ohio State University. This code, called the Antenn a Workbench 



(iKouvoumiian & Path^[l974 



Lee etal.ll 1979: 



Lee & Ruddick 1985 



Lee et al.lll99C ), can be used 



to analyze single or multiple reflectors. It is capable of calculating both near-field and far-field 
radiation patterns. The theoretical approach is based on a combination of Geometrical Theory of 
Diffraction (GTD) and Aperture Integration (AI) techniques. Typically, AI is used for computing 
the main beam and the near sidelobes while GTD is used for computing the far sidelobes includ- 
ing the backlobes. For near-field calculations, GTD is used in some cases for the whole region 
including the near axis region. In addition, the code has the conventional Physical Optics (PO) 
method option, where the currents on the reflector surface over a two-dimensional grid system are 
integrated to obtain the radiation pattern. The code usually uses the PO option for antennas with 
multiple reflectors. The Antenna Workbench was modified specifically for the large main reflector 
of the GBT. In the present analysis, the multiple reflector option using the PO technique was not 
used because it does not include edge effects and hence cannot compute the far sidelobes accurately. 

Fig. lA.ll shows the geometry of the primary and secondary reflectors of the GBT in the antenna 
symmetric plane that is vertical given the Alt-Az mounting. Three coordinate systems to be de- 
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Fig. A.2. Subreflector near-field patterns in the "s" coordinate system as a function of angle ds from 
the ellipsoid axis Zs, in half-planes with different values of (p/. (ps = -90°, red; 0°, green; 90°, blue. 
Vertical lines indicate the near edge (9°.4), center (54?0), and far edge (87°4) of the main reflector 
in flie ^ = -90° half-plane. 



scribed further below are indicated. The non-subscripted one is associated with the primary, for the 
aU-sky power pattern. That associated with the subreflector is subscripted A third coordinate 
system subscripted "p" is used for subreflector illumination of the primary, including for calcula- 
tions of spillover past the primary. These are right-handed x-y-z systems, with the x axis directed 
up out of the page and yz in the antenna symmetric plane. In these coordinate systems we define 
(p to be the azimuthal angle (range -90° to 90°) around the z axis, measured relative to the x axis 
((f> = 0°) toward the y axis ((p = 90°). The polar angle 9 is measured in half-planes of constant (p 
with an edge along z, relative to the z axis {0 = 0°) and ranging over 180° to the -z axis. 

The radiation pattern of the GBT was computed in two steps. First, the pattern of the subreflec- 
tor was calculated, in the coordinate system whose origin is at one of the foci of the ellipsoid 
of which the subreflector is a segment. This focus is also the primary focus (PF) of the main reflec- 
tor (Fig. lA. Il l and is the phase reference point for the subreflector patterns. The L-band feed horn 
is at the other focus of the ellipsoid, which is the GBT Gregorian (secondary) focus (SF). The Zs 
axis runs f rom PF toward SF. These calculations used the measured far-field patterns of the L-band 



feed horn (iSrikanth 



19931) to illuminate the subreflector. The grids on the subreflector were set at 



5x5 cm. Since the main reflector is in the near field of the subreflector, the near-field option in the 
code was used. 

Figure lASl shows the power patterns calculated for the subreflector in five 0, half-planes. Given 
the orientation of the "i" coordinate system, it is the patterns in the half -planes approaching (ps - 
-90° that illuminate the main reflector. In this coordinate system, in the antenna symmetry plane 
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Fig. A.3. Subreflector illumination patterns on the main reflector as a function of angle 9p from Zp, 
the direction from PF to the center of the primary, for three values of (pp in the main reflector ("p") 
coordinate system: (pp - -90°, red; 0°, green; 90°, blue. Vertical lines indicate the near (44°6) and 
far edges (33?4) of the main reflector in the - +90° half-planes, respectively. 

the near edge of the main reflector is at 6s = 9°4 and the far edge is at 87?4. At 54° is the center 
of the main reflector]^ The power pattern shown for cf>s - -90° is appropriately peaked to provide 
good illumination. 

To calculate the illumination of the main reflector and the spillover past its edge, these pat- 
terns were transformed from the "s" to the coordinate system with origin at PF and axis Zp 
oriented from PF to the projected center of the 100-m main reflector ("Fig. lA.lt . Fig. lA.3l shows the 
illumination of the main reflector from the subreflector transformed to this coordinate system, for 
three diflferent (pp half-planes. The edge of the main reflector subtends an angle that varies between 
Op = 33°4 and 44°6 as (pp changes from -90° (far edge) to 90° (near edge). As can be seen from 
both Figs. lA.2] and lA.3l the illumination taper in the antenna symmetry plane at the far edge of the 
main reflector is about 5 dBi and somewhat lower, -13 dBi, at the near edge, so that the spillover 
lobe from the main reflector is not symmetrical around Zp (see below). 

The second step was to determine the all-sky response of the GBT, expressed in the 
non-subscripted coordinate system in Fig. I A. II Subreflector patterns as calculated by Antenna 
Workbench for 72 half-planes in the range of -90° < (p,, < 90° were used. The main reflector 
was gridded into regions 30 x 30 cm, which is about 1 AA, adequate for the present work using AI. 
Far-field patterns were calculated at (f.Q5 intervals in <p and 0°02 intervals in 6. The computation 

This is the projected center as seen along the direction of the main beam, z; the projected aperture diameter 
is 100 m. 
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Fig. A.4. Calculated all-sky response of the GBT at 1 .4 GHz, the main beam and far sidelobes, as 
a function of angle from the main beam, 6, for three values of (p in this coordinate system: (p - 
-90° (towards the feed arm in the antenna symmetry plane), red; 0° (perpendicular to the antenna 
symmetry plane), green; 90° (away from the feed arm), blue. To show details of the sidelobes the 
scale has been magnified so that the forward gain of the main beam, 61.6 dBi at = 0° (Fig.|5] 
Table[TJ, is beyond the range of this figure. Vertical lines indicate the far edge (98?2), center (13 1?5), 
and near edge (176?2) of the main reflector in the (p - 90° half -plane. 

of this pattern switches from AI to GTD 6 - 2°65 = sin"'(l/ '</A^) where A„, is the aperture 
diameter in wavelengths in any given (p plane. 

Fig. lA.4l shows the main beam and the far sidelobes of the GBT calculated for three diff'erent (p 
half-planes. The on-axis (forward) gain is 61.61 dBi. An expanded view of the response near the 
main beam is given in Fig. |5] The calculated half-power beam width at 1.40 GHz is 9! 17 x 9! 14, 
which differs from the measured values by < 1% and < 2%, respectively. 

The sidelobes stay mostly below -20 dBi. The peaks at about 10 dBi seen for (p = 90° at 5° 
and for = 0° at 13° are the spillover sidelobe from the subreflector (Appendix |BT3]). This can be 
appreciated in Figs. H] and |9l noting that the -V axis there corresponds to (p - 90° and the H axis 
to = 0°. 

The (p = 90° plane is in the symmetric plane of the GBT structure in the direction away from the 
feed arm and so as increases the main reflector is crossed. The two sidelobe peaks near 96° and 
178° are caused by the spillover of the subreflector illumination (Fig. lA.3l l past the far edge (98°2) 
and near edge (176°2) of the main reflector, respectively. Note that the former is stronger because 
as shown above the illumination from the subreflector is larger at the far edge. This spillover lobe 
is manifested as an asymmetrical ring around Zp as (pp changes. The peak near - 131°, which 
lies close to Zp (131°5), is caused by the diffracted rays off the main reflector edge, similar to the 
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backlobe that is observed on an on-axis antenna. This is like the Arago spot from the subreflector 
discussed in Appendix IB. 21 Note that it is offset from the cone axis defining the main reflector 
(137?2) to lower 0, toward the far edge, because of the tilt of the reflector as seen from PF. 



Appendix B: Calculation of diffraction from the subreflector using Fresnel 
diffraction theory 

The primary reflector of the GBT is an offset seg ment of a symmetric parabola designed to elim- 
inate aperture blockage and attendant sidelobes (IPrestage et alJl2009l) . The Gregorian secondary 
subreflector is an offset section of an ellipsoid with the prime and secondary foci of the GBT at 
its two foci (Fig. lA.lb . The GBT L-band receiver feed horn is at the secondary focus. The geo- 
metric path from feed to subreflector to main reflector is clear of any obstruction. However, to a 
receiver feed that does not have a perfectly sharp cutoff of its field pattern at the subreflector edge, 
the subreflector appears as a blockage of the sky and so its diffractive eff'ects must be considered. 
The GBT L-band feed does indeed somewhat over-illuminate ("spills over the edge of") the sub- 
reflector, because a feed having a more sharp ly tapered beam pattern would be too large and heavy 
for the feed turret ( Norrod & Srikantlilll996 ). Therefore the field pattern of the feed couples to the 
diffracted radiation, producing a modified Fresnel diffraction pattern, whose main "edge" feature 
is a set of annular peaks (rings) centered on the subreffector and beyond its projected edge. The 
dominant first ring is often called the "spillover sidelobe." Here we present a calculation of the full 
power pattern. 



B.1. Simple diffraction pattern of ttie subreflector 

The shape of the subreffector is defined by the intersection of a cone, with apex at the secondary 
focus and an opening half-angle 6h - 14.99°, with an offset ellipsoid of revolution (a prolate 
spheroid). As seen from the feed, the subreffector is to a first approximation a circular blockage 
defined by radius a - 3.775 m at distance d - 14.1 m from the secondary focus. At the oper- 
ating wavelength of 21 cm, the secondary focus is in the near-field (d < a^/4A) with respect to 
the subreffector and therefore satisfies the Fresnel condition. Furthermore, a" jd^^A < 1 so that 
the usual quadratic expansion of the phase used in the "Fresnel approximation" is valid (e.g.. 



Sheppard & Hrvnevvch 



1992 : 



Daugerlll996l) . Thus the characteristics of the sidelobe pattern are 



closely related to the familiar near-field Fresnel diffraction by a circular blockage. This can be cal- 
culated for a receiving system with a plane wave arriving along the cone (negative Zc) axis or, by 
reciprocity, for a transmitting system with an isotr opic source at th e feed. We have evaluated this 
from the analytic solution for a receiving system by 
numerically using the phasor approach described by 



Hovenad(ll989h and obtained the same results 

icit 



Daugei ( 19961) . in which reciprocity is impl 
in the calculation. Dauger's Fresnel Diffraction Explorer (FDeJ can also be used to explore this 
interactively. 

A radial profile of the axially-symmetric power pattern is shown in Fig. IB. II In this pattern, and 
in the GBT sidelobe pattern, there are two main phenomena of interest: (i) the finely-spaced ring 
structure near the cone axis with related modulations continuing at angles beyond the edge of the 



http://daugerresearch.com/fresnel 
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Fig. B.l. Radial profile of the familiar diffraction power pattern of a circular blockage, 6 being with 
respect to the axis of the cone defining the subreflector. The vertical dashed line marks the edge 
of the subreflector as seen from the secondary focus. The dashed curve appro ximating the pat tern 
near the central Arago spot is the lowest-order term in the analytic solution bv iHovenad d 1 9891) . 



subreflector, and (ii) the coarser ring structure existing only beyond the edge. The existence of the 
diffraction spot on axis depends on interference of waves arising from all around the circumference 
of the subreflector in this circularly-symmetric geometry. On the other hand, the outer coarse ring 
structure is an "edge effect," depending only locally on the sharp boundary, and so is qualitatively 
similar for a circular, square, or knife-edge configuration. As discussed below, the actual GBT 
sidelobe pattern shows the effects of three other factors: (i) the tapered beam of the L-band feed at 
the secondary focus, (ii) a screen (Fig.O adjacent to the subreflector edge in front of the supporting 
feed arm to redirect feed spillover, and (iii) the tilt of the subreflector. 



B.2. The Arago spot 



In the center of the diffraction pattern, at - 0, is the Arago (or Poisson) spot, surrounded by a 
series of fine rings of c omparable width . This pattern is described accurately by the lowest-order 
term in the solution bv iHovenad (1 1989b . i.e., the square of the zero-order Bessel function of the 
first kind, J^^{v), where v = kar/d, k - In /A, r is the radial coordinate in the image plane and 
- arctan(r/ii). For the inner rings, the positions of the minima and maxima are near the zeros 
and peaks of the function 1 + sin(2v), and thus at (4« + zp)in IA)lka, n - 1, 2, ... and zp - -1 or 1, 
respectively. The first maximum is at \°.99. In these units, the half-power point of the main peak, at 
V = 1.1263, is at angle 1.434(:7r/4)/^a,andsotheFWHMof the Arago spot is 1° 14; this is inversely 
proportional to a and, for fixed cone angle, also d. The predicted FWHM and ring positions are in 
close agreement with the properties of the actual scans of the GBT sidelobe using the Sun. 
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Fig.B.2. Solid (black): simulation of the effect of the tapered beam on the diffraction pattern, 
normalized to unity at the Arago spot. Step function (red): tapered field pattern of the L-band feed, 
relative to unity at the edge of the subreflector. Dashed (green): radial profile in the H direction 
through the observed sidelobe pattern. 



B.3. The main spillover sidelobe 



The dashed curve in Fig. IB. 21 shows the actual radial profile of the observed sidelobe through the 
Arago spot along the H direction. The oscillating coarse ring structure beyond the subreflector is 
an edge effect. However, the angular extent and amplitude are not the same as the pattern shown 
in Fig. IB. II because the strongly tapered feed pattern, already down 14.7 dB at the edge of the 
subreflector, does not illuminate the region beyond the edge of the subreflector uniformly. Thus the 
contribution to the diffraction pattern from radiation at large angles is strongly suppressed. We have 
simulated this as a steadily increasing "blockage" beyond the ed ge of the subreflector by imple- 



menting a method for treating variable transmission suggested bv lDaugerl (119961) . The amplitude 
and phase of the far-field pattern have been measured for the prototype L-band feed ( ISrikanthll993h 
and as assumed in Appendix lAl should be approximately valid for the distance of the subreflector. 
The phase is roughly constant and we have taken the transmission (coupling) to decrease following 
the field pattern. The resulting diffraction pattern is shown in Fig. IB. 21 normalized to the amplitude 
of the Arago spot. The outer sidelobe pattern is systematically and differentially diminished: the 
peak intensity of the major spillover lobe compared to the Arago spot is lower and the outer rings 
are suppressed even more. As expected, the position of the peak of the major ring moves closer to 
the edge of the subreflector. These effects are evident in the observed sidelobe pattern, although the 
details are somewhat different; this might arise in part because of the large angles encountered in 
the actual GBT configuration. 

B.4. Effects of the screen 

The diffraction pattern is altered by the screen (Fig. [3]l because it eliminates the contributions 
of interfering waves from a wedge with angular range <ps ~ 40° along the circumference of the 
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subreflector centered on the positive V axis (away from the main beam). The amplitude of the 
Arago spot is reduced by a factor (1. - (ps/36Q)^ ~ 0.8. In the simulation for Fig. IB. 21 we have 
not included the screen; this would lower the Arago spot while leaving the major spillover lobe 
unaffected along the H profile thus bringing the relative deflections even closer to that observed. 

Our calculations incorporating the screen indicate that in the profile along H the first fine ring is 
enhanced but in the orthogonal direction, along V, the ring amplitudes are reduced while the radial 
offset of the rings increases. These more subtle effects are discernable in the observations. 

Because the spillover sidelobe is an edge effect, the disruption of the edge of the subreflector 
by the screen has a profound effect locally. This is clear in the observed pattern (see Figs.[8]and|9]l 
where there is a bite out of the spillover ring pattern, about a 40° wedge out from the Arago spot 
at positive V centered on the feed leg. The screen apparently acts as an edge too, brightening the 
adjacent portion of the spillover sidelobe. 

Superimposed on the inner edge of the observed main spillover lobe at negative V near H - 0, 
on the side of the main beam opposite to the direction of the screen, there is a pair of features 
in excess of the simple diffraction pattern (see Figs.[8]and[T0t. These have about the same beam 
integral as the portion of the major spillover lobe removed in the wedge, qualitatively consistent 
with conservation of energy. The screen has been well designed so that these features too lie well 
away from the main beam. 



B.5. Effect of the tilt and motion of tine subrefiector 



The required GBT geometry dNorrod & Srik:anthlll996l) is such that the subreflector is not quite 
perpendicular to the axis of the cone, but rather tilted by t = 17?0; relative to the defining cone axis, 
the most distant edge is toward the main reflector and the direction of the main beam (Fig. lA.ll ). 
The semi-major axis, in the antenna symmetric plane and joining the edges of the subreflector along 
the direction of tilt, is of dimension 3.975 m, and the semi-minor axis is a = 3.775 m. The distance 
along the cone axis from the secondary focus to the center of the minor axis is then d - 14.1 m. 
These were the values of a and d used above. The tilt shifts the position of the Arago peak relative 
to the direction to the subreflector center, such that there is still constructi ve interference between 
waves from the distant and near edges. Solving Eq. A7 of iDaugeii ( Il996h . phases for the ends of 
the major axis are equalized for an offset 2?05 tan t in the direction toward the main beam. For the 
specified tilt, our estimate of the offset is (f.63. The angle between the nomin al cone axis defining 
the subreflector center and the main beam is 12?329 JNorrod&Srikanthll996i at the rigging angle, 
elevation 50°.3 ( Nikolic et alj2007 ), the geometry used in the above calculations. Thus the expected 
offset of the Arago spot at the rigging angle would be 1 1°.7 along the V axis above the main beam. 

At the rigging angle the main reflector is a perfect paraboloid to 21 -cm radiation and the feed 
and subreflector are located optimally. At any other elevation angle, the main reflector deforms, 
displaces, and rotates while the feed and subre flector displace later ally in the symmetry plane, 
resulting in a slight misalignment of the optics jSrikanth et al.lll994l) . To track the primary focus 
position and maintain optimal efficiency together with precision pointing, the subreflector is moved 
to compensate (see, e.g., GBT Commissioning Memos 7 and ijfl). Thus the angle between the 



http://www.gb.nrao.edu/ rmaddale/GBT/Commissioning/memolist.html 
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cone axis of the subreflector and the main beam changes sHghtly while tracking a source during an 
observation. 

The elevation of the telescope when the Arago spot crossed the Sun during two different nearly 
vertical scans was 3 1°5 and 35°.4. Allowing for the change with elevation we estimated that the 
position of the Arago spot would be at 1 1?8. The two scans gave identical observed offsets, 1 1°.9. 
The good agreement is perhaps fortuitous, given the complexity of the calculation. 

B.6. Comments on the derived sidelobe model 

These simple simulations reveal potential subtle issues regarding our method of deriving the side- 
lobe pattern by scanning it with the Sun. One involves the data reduction, where of necessity the 
scans have had a baseline removed. This introduces some uncertainty in the derived pattern. For ex- 
ample, in the area of the Arago spot there is a plateau not explained by the above simple diffraction 
theory (Fig. lB.2l i. Even though this plateau is relatively weak, its subtended area could potentially 
lead to a significant contribution to the stray radiation. Another possibility, with an effect of the 
opposite sign, is that the derived pattern does not sink quite to zero in between the main spillover 
lobe and next minor rings. Again, because of the large area, even weak perturbations like this could 
have considerable power. To compensate for inadequate baseline removal, an approach could be to 
introduce a smooth parameterized sidelobe contribution centered on the subreflector and a function 
only of 9, e.g., a low order polynomial plus exponential. The parameters could be optimized during 
the procedure described in § l6.1.2l to determine the optimal scaling of the sidelobe pattern. The ef- 
fective coverage of different parts of the sky by the sidelobe pattern would therefore be reweighted. 

A second issue concerns how the actual diffraction pattern might change in position with ele- 
vation of the telescope away from the rigging angle. As discussed in connection with the position 
of the Arago spot, displacement of the cone axis of the subreflector with respect to the main beam 
is difficult to calculate precisely but it might be several tenths of a degree over the range of eleva- 
tions used in our survey, from 15° almost to the zenith. A smaller range was used in our mapping 
of the sidelobe pattern with the Sun. Because major features of the sidelobe pattern are related to 
subreflector spillover, and the latter is geometrically tied to the cone axis, the whole pattern can 
shift slightly with elevation. As a further complication, the screen is fixed to the feed arm, not 
moving with the subreflector during focus tracking. Thus the sidelobe pattern as measured and im- 
plemented will not be perfect; this might account (at least in part) for the slight differences among 
the optimal values of ?7si derived independently in different survey regions (Sect. 16.1.31 ) and for the 
stray radiation correction errors — which are nevertheless small (Sects. 1731 and l7.4] ). 

Appendix C: Stray radiation correction algorithm 

The stray radiation correction is computed by evaluating the integral in Eq. |4l This is based on 
a program written in 1997 by Edward M. Murphy to correct stray radiation for HI spectra from 
the 43 m (140 ft) telescope at the NRAO incorporating a multi-component semi-analytic model 
developed from careful measurements of the 43 m telescope and its sidelobes (E.M. Murphy, 1993, 
unpublished M.A. thesis. University of Virginia). For each input H I spectrum, an output spectrum 
of the estimated stray radiation was produced, with the input file serving as a template for the 
output file. The algorithm to estimate the sidelobe was straightforward, using an empirical model 
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of the H I sky augmenting an early version of the Leiden-Dwingeloo survey with data from the Bell 



Labs survey and the Parkes 60-foot survey jKerr et al.ll 1986 ; 



Stark etalJI 1992 : 



Hartmann & Burton 



19971). This model sky had HI spectra at half-degree intervals in Galactic latitude / and longitude b. 



For each (/, b) point in the model, the date and time of the input spectrum were used to determine the 
corresponding position on the sky as viewed from the telescope. If this point lay above the horizon 
(it is assumed that there is no HI emission reflected from the ground; see Appendix lD.3l ) and was 
more than one degree from the main beam, then the corresponding survey spectrum was multiplied 
by the relevant solid angle, by a factor to account for atmospheric absorption at the corresponding 
zenith angle, and by the sidelobe level at that position relative to the main beam. The appropriate 
velocity shift relative to the LSR was applied and the spectrum was spline-interpolated onto the 
input-spectrum velocity grid. The sum of all such weighted velocity-shifted spectra, looping over 
the grid of the model H I sky, was output as the stray radiation estimate. This computation of one 
stray-radiation spectrum required about one minute of CPU-time on a typical workstation. 



C.1. Modified algoritiim for the GBT 

Major modifications to support calculation of stray radiation corrections specifically for the GBT 
involved four main areas: adding support for the SDFITS file format used for GBT observations, 
modifying the sky position algorithms to account for the fact that the GBT has an altitude-azimuth 
mount instead of the equatorial mount of the 43 m telescope (as well as including the horizon as 
seen from the GBT), modifications to speed up the calculation by a couple of orders of magnitude, 
and finally obtaining a good map of the GBT sidelobes (Sect.|6]l which is stored in and interpolated 
from a 2-D array. 

For the empirical model H I sky we adopted the LAB all-sky survey (iKalberla et al after 



first subtracting a Gaussian with peak amplitude 0.048 K, center at -2 2 km s ' , and FWHM - 
167 kms ' from each LAB spectrum (the so-called "Wakker correction". IWakker et al I2OI1I) and 
then dividing by a scale factor 1.0288, as discussed in Sect. 18.21 (see Eq.|9]l. The "Wakker correc- 
tion" to the LAB survey has a very small effect on the calibrated GBT spectra as it enters scaled by 
T/si times the fraction of rj^x that is above the horizon. As the product of these quantities is always 
< 0.1, its effect on calibrated GBT spectra is always < 0.005 K. 

Air mass measurements were available for NRAO at 5° intervals in elevation angle el for 0° < 
el < 90° and at 1^5 intervals for Q° < el < 15° (R.D. Maddalena 2005, private communication: 
the air mass measurements were actually performed nearby at Hot Springs, VA, on 15 June 2005). 
These measurements were used to obtain the atmospheric optical depth Tatm as a function of el: a 
zenith optical depth at 1420 MHz of Tzenith ~ 0.01, and an air mass of roughly 31 at the horizon, 
implies that Tatm ~ 0.32 at el - 0° (rather than becoming infinite there). However, as the sky at 
el < 1 ° is obscured by hills around most of the horizon, the improved atmospheric absorption values 
had little effect on the calculated stray radiation. Information on the exact weather conditions is not 
used at present when calculating the opacity, though there are indications in our data that rain can 
add an extra opacity that might be accounted for in the future (Sect. l7T4l R.D. Maddalena, private 
communication). Snow might also have a significant effect both with its precipitation and with its 
accumulation in the GBT dish — any effects from the latter not being sufficiently predictable to be 
corrected for. However, variations in atmospheric pressure, temperature, and humidity are expected 



46 



Boothroyd et al. : Accurate GBT H I Measurements 

to have very little effect on the atmospheric absorption; this is confirmed by the fact that the repeated 
observations of S6 and S8 (Sect. l7!4l) show no evidence for a systematic summer- vs. -winter effect. 

Refraction results in more of the sky being visible to the telescope, an extra strip approximately 
(f.3 wide around the horizon. Effects of temperature, pressure, and humidity are included in the 
refraction calculations. Inclusion of refraction typically changes the calculated stray radiation by 
less than 1%, but up to 10% in extreme cases. 

For integration times of a few minutes or more, there are significant changes in the position of 
the sky relative to the horizon and, for the GBT, in the rotation of the sidelobe beam relative to the 
sky with accompanying changes in the Doppler shift of the stray radiation being accumulated. For 
an integration time of 120 sec, such as we used for "staring observations" as opposed to on-the-fly 
mapping, calculating the sky, horizon, and beam positions at the middle of the time interval (rather 
than the beginning) yields typical improvements of a fraction of a percent though in extreme cases it 
can be 10%. To allow for even longer integration times, calculation of stray radiation contributions 
at multiple times within a single observation is supported. 

The calculation was sped up mainly by two modifications. First, speed-up by nearly an order 
of magnitude was obtained simply by using linear rather than spline interpolation in velocity for 
the spectrum being accumulated. Second, another order of magnitude was obtained by introducing 
a four-level tiling tree for the model sky spectra: a four-level tesselation of the sky. The deepest 
level corresponds to independent spectra in roughly (f.5 x 0°5 areas on the sky, appropriate to the 
LAB survey sampling. Near the Galactic equator, these tiles are at 0^5 intervals in b and /, i.e., 
the LAB grid. Near the Galactic poles, LAB spectra at multiple / values were averaged together 
to create tiles with widths on the sky of at least 0°25 but no more than 0°5. The next level up, 
coarser tiles of roughly 1° x 1°, comprised of averages (weighted by solid angle) of four (or, near 
the poles, sometimes three) of the lowest-level spectra. Similarly, two further higher tiling levels 
of roughly 2° x 2° and 4° x 4° were obtained. The loop over the tiles of the model starts at the 
highest 4° x 4° level. Since much of the sidelobe area is quite smooth, this tiling made it possible 
to reduce greatly the number of spectra that had to be accumulated for each stray calculation. 
Where the sidelobe varied significantly on length scales smaller than the relevant tile, or where 
part of the tile would be below the horizon or adjacent to the main beam, the appropriate tiles on 
the levels below were traversed instead. The magnitude of the introduced errors was estimated by 
computing stray radiation corrections for a random sample of positions, dates, and LST values. 
For these modifications that yield faster sidelobe spectrum evaluation, the resulting errors are small 
(rms 0.2%, max 3% of the stray correction). 



Appendix D: Uncertainties not related to the GBT response pattern 

D.1. Undersampling in the LAB survey 

The correction for stray radiation relies on the LAB survey to provide a 21 cm T\, over the entire 
sky visible from Green Bank. The main contribution to the LAB survey for the part of the sky 
seen by the GBT is the Leiden-Dwingeloo (LP ) HI survey that consists o f observations spaced 



(K5 on the sky at an angular resolution of 35' (iHartmann & Burton 



1997). It is sometimes said 



that the undersampling makes the LD survey equivalent to a survey with an angular resolution 
of order 1°, but this is not correct. Because of the undersampling, the LD survey (and hence the 
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LAB in the North) cannot be used for accurate interpolation between the measured positions - 
information at spatial frequencies that are resolvable to the Dwingeloo antenna are aliased into 
lower spatial frequencies biasing our estimate of the H I sky. At present there is no remedy for 
this and it contributes an unknown error to the stray radiation correction. A new H I survey of the 
northern sky is underway that should provide a m uch more accurat e, complete, data set from which 



to calculate the stray component in GBT spectra (IKerp et al.ll201 II) . 



D.2. Errors in the LAB survey 

Any errors in the LAB survey spectra will propagate through the stray radiation correction, reduced 
in their effect by t/si = 0.098 1 and the fractional solid angle of the erroneous spectra, and by the fact 
that multiple errors in op posite directions m ight cancel out. We consider here errors in the LAB 
survey discovered earlier (iHiggs et al and during the course of this work. 



As discussed in Sect. 18.21 a very few of the LAB survey spectra have normalization errors of 
up to about 30%, but only a small fraction of the spectra appear to have normalization errors even 
at the ten percent level. 



As noted earlier bv lHiggs et al.l (120051) when they compared their DRAG 26-m data to the LAB 
survey, a very few LAB spectra have large spurious features. These are detected by comparison 
of LAB with other data; the fact that it is the LAB feature that is spurious is indicated by the 
fact that (1) the feature appears in the LAB spectrum but not in the other data, (2) nor does the 
feature appear in LAB spectra at adjacent latitudes, and (3) the feature often has a shape that 
is not typical of any real feature. (Note that LAB survey spectra at adjacent longitudes are not 
independent at high latitudes, where a spurious feature might thus propagate with varying intensity 
across several longitude points.) Gnly where alternative data are available can the above errors be 
reliably detected and corrected. This alternative is available for only 13165 of the 257762 LAB 
spectra; in 125 of these 1 3165 spectra, a tot al of 141 spurious features were found (including 33 



in the regions surveyed bv lHiggs et al.ll2005h . Most were small, "dips" or "peaks" with a width of 
about 2 km s ' superimposed on real features an order of magnitude larger in both width and height. 
However, 6 of the LAB spectra had very wide spurious "peaks" at negative velocities. The LAB 
spectrum at (I - l3K'.5,b = 37°) has a "peak" over the range -300 kms"' < v < -160kms~^ 
with a plateau of height 3.6 K in the range -240 km s ' < v < -190kms ' (this region should 
be zero); this feature also shows up at / = 139° and 139?5, albeit with lower amplitudes of 0.5 
and 0.1 K, respectively. The LAB spectrum at (/ = 141°,/? = 41';5) has a peak of height 0.6 K 
in -190kms ' < v < -70kms ' (this region should be zero). Finally, the two LAB spectra at 
(I = 132° and 132°5,/7 = 47°) have a peak of height 2 K in -240 kms"' < v < -20kms-i 
(the convolved GBT spectra, and adjacent LAB spectra, have three narrower peaks in this velocity 
range, with heights of about 0.2 K). Replacing the spurious features by a linear interpolation over 
the relevant velocity range works very well, except in the last two cases (where convolved GBT 
spectrum values can be used instead, in the relevant velocity region). 

The existence of spurious LAB features and random LAB calibration errors in our large GBT 
fields allows for some quantification of the potential propagated error. For pointings near the middle 
of these fields, the major contribution to the stray radiation comes from the LAB spectra that were 
corrected for the above two types of errors, and comparing these cases to computations without 
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these corrections indicates that these spurious features and mis-calibrations introduce errors of at 
most a few percent in the computed stray radiation correction. Thus any such spurious features 
that cannot be corrected (i.e., in regions with no comparison data) are not expected to contribute 
significant errors to the stray radiation computation. 



D.3. Errors from the stray-correction algorithm 

The program currently assumes that no 21 -cm radiation is picked up through reflection from the 
ground; the tree-covered mountains around the GBT can be assumed to be good absorbers, but 
the fields in the immediate vicinity of the telescope might have a non-t rivial albedo especiall y 
when wet. This was found to be the case for the moors around Dwingeloo ( Hartmann et alJ[l"996 ). 
However, as the dominant GBT sidelobes are near and above the main beam (Fig. [8]l, the fraction 
of the response lying on the ground is at most a few percent and so the effect of reflections should 
be quite smaU. 

The "sidelobes" over which the stray radiation is accumulated are considered to begin a degree 
from the center of the beam. Shifting this cutoff from a 1.0° radius to 0.9° has an effect of rms 2% 
(max 14%) on the stray radiation, according to a set of random test cases. Such effects are expected 
to be most significant when the GBT beam is pointed at a region with extended strong signals, 
such as the Galactic plane. For pointings in regions of low signal, even much larger changes in the 
definition of where the sidelobes begin have little effect. 

As discussed above (Sect. lB.6T l. variations in the geometry of the optics with telescope elevation 
angle can shift the sidelobe pattern by several tenths of a degree (its shape and amplitude pattern 
might also change slightly). This is at least an order of magnitude smaller than the angular scale 
on which the sidelobes vary, but could nonetheless lead to significant stray radiation errors for 
cases where the mean sky brightness is varying significantly at the edge(s) of the spillover lobe. 
Typically, such errors would be expected to be no more than a few percent of the calculated stray 
radiation; but in extreme cases, this effect could be larger. The present stray-correction algorithm 
ignores this telescope elevation effect, but in principal measurements might allow at least a partial 
correction for it. 

For pointings in bright complex regions, spectra obtained at different times might change by a 
few percent due to the rotation relative to the sky of the slightly asymmetric main beam and inner 
sidelobes within 1.0° radius. Such "near-beam" effects cannot be corrected by this program due to 
the coarse grid of the all-sky LAB survey. This could be approached as a deconvolution problem 
using the accurately known near beam if a small fully sampled map were obtained quickly with the 
GBT. Given our interest in high latitude emission, this has not been pursued here. 
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